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As a Direct-Reading Instrument The extreme sensitivity 
of the d-c amplifier is utilized to check plasticizer insula- 
tion resistance values in the megamegohm range. 


As a Recorder Preamplifier The rack-mounted ampli- 
fier above is being used to increase the sensitivity of a 
recorder in running special tests of switches. 


As a Null Detector The d-c Amplifier is being used above 
for factory checking and calibration of instruments. 


L&N Low Level 
d-c Amplifiers are 


@ In response (o the constant demand for versatility 
in precision instruments, these d-c Indicating Ampli- 
fiers combine the functions of three useful instru- 
ments in one: 


1. A Direct-Reading Instrument that is al- 
ways ready to use... never any readjusting of 
zero, either initially or ‘during a series of 
Simply select the range in which you want to 
work by turning scale-multiplier knob. 


2. A Recorder Preamplifier — Values meas- 
ured by Stabilized d-c Amplifiers can be recorded 
directly on Speedomax recorders. 


3. A Null Detector more sensitive than most 
ters, yet with full scale re- 
sponse time of only 2 to 3 seconds. These instru- 
ments are unaffected by vibration; leveling is 
not necessary. At the turn of a range knob, a 
wide choice of sensitivities can be obtained with- 
out external shunts. A non-linear response char- 
acteristic is also available for easy balancing. 


These amplifiers are suitable for handling low level 
measurements with thermocouples, strain gages, 
bolometers—bridge and potentiometer circuits— 
ionization, leakage and phototube currents—almost 
any measurement of extremely small direct current 
or voltage. 

Self-contained, the unit can either be used “‘as is” 
or removed from case and mounted on a 19” relay 
rack. 

For details, including complete specifications, send 
for Folder EM9-51(1). Write our nearest office, or 
4992 Stenton Ave., Phila. 44, Penna. 
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REMOTE LIQUID. 
LEVEL INDICATOR 


Now ...a more accurate check on hard-to- 
see overhead boiler water gages... right in 
front of your eyes on the instrument panel or 
other convenient place. 

Operated by the boiler water itself... by the 
pressure differential between a constant head 
of water and the varying head of water in 
the boiler drum. 

Indicating mechanism is never under pres- 
sure, due to unique magnetic transmission 
from pressure side to indicating side of instru- 
ment. Action is instant, constant, frictionless. 
There are no stuffing boxes. Mechanism is per- 
fectly balanced on jewelled bearings outside 
of the pressure chamber. Suitable for all 
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Deep in the ocean, a submarine 
telephone cable system is extremely 
hard to get at for adjustment or re- 
pair. This makes it vitally important 
to find out what can happen to such 
a system before it is installed. 


Bell Laboratories engineers do this 
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ocean floor conditions on dry land. 
Among many factors they test for are 
the effects of immense pressures on 
amplifier housings and their water- 
resistant seals. They also test for 
agents which work very slowly, yet 
can cause serious destruction over the 
years—chemical action, marine borers 
and species of bacteria which strange- 
ly thrive under great pressures. 


Through this and other work, Bell 
Laboratories engineers are learning 
how to create better deep-sea tele- 
phone systems to connect America to 
the rest of the world. 


Highly precise instruments developed by Bell Lab- 
oratories engineers are used to detect infinitesimal 
changes in cable loss—to an accuracy of ten mil- 
lionths of a decibel. 


“Dry Land Ocean,” under construction at Bell Labora- 
tories, simulates ocean floor conditions, is used to test 
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pleted trough is roofed in and is filled with water which 


maintains the pipes at 37° F., the ocean floor temperature. : 
Seawater and sediment in bottle characterize 


ocean floor. Test sample of insulation on coiled 
wire is checked for bacterial attack by conduct- 
ance and capacitance tests. 
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inventions of high order and for particularly meritorious improvements and developments in 
machines and mechanical processes. In the event of an accumulation of the fund for medals 
beyond the sum of one hundred dollars, it is competent for the Committee on Science and the 
Arts to offer from such surplus a money premium for some special work on any mechanical 
or scientific subject that is considered of sufficient importance. 

The Louis E. Levy Medal (1923—Gold Medal).—This medal is awarded to the author 
of a paper of especial merit, published in the JouRNAL OF THE FRANKLIN INSTITUTE, prefer- 
ence being given to one describing the author’s experimental and theoretical researches in a 
subject of fundamental importance. 

The George R. Henderson Medal (1924—Gold Medal).—This medal is awarded for 
meritorious inventions or discoveries in the field of Railway Engineering. 

The Walton Clark Medal (1926—Gold Medal).—This medal is awarded to the “author 
of the most notable advance in knowledge or improvement in apparatus, or in method con- 
cerning the science or the art of gas manufacture or distribution or utilization in the produc- 
tion of illumination, or of heat, or of power.” 

The Frank P. Brown Medal (1938—Silver Medal).—This medal is awarded to in- 
ventors for discoveries and inventions involving meritorious improvements in the building and 
allied industries. 

The Newcomen Medal (1943—Gold Medal).—This medal is awarded, not oftener than 
once in three years, for achievement in the field of Steam. 

The Francis J. Clamer Medal (1943—Silver Medal).—This medal is awarded at least 
ence in five years for meritorious achievement in the field of Metallurgy. 

The Stuart Ballantine Medal (1946—Gold Medal).—This medal is to be awarded in 
recognition of outstanding achievement in the fields of Communication and Reconnaissance 
which employ electromagnetic radiation. 

The Boyden Premium (1859).—This premium is awarded not oftener than once in five 
years to any resident of North America who has recently made a notable experimental deter- 
mination of the speed, in free space, of radiation in any region of the entire spectrum. 


The William M. Vermilye Medal (1937—Bronze Medal).—This medal is awarded not 
oftener than biennially in recognition of outstanding contribution in the field of Industrial 
Management. 


The Certificate of Merit (1882)—A Certificate of Merit is awarded to persons ad- 
judged worthy thereof for meritorious inventions, discoveries or improvements in physical 
processes or devices. 

For further information relating to these awards apply to the Secretary, Committee on Science and the Arts 


- 

2 

: 

id 

: 

= 

; ae 


JourRNAL OF THE FRANKLIN INSTITUTE 


Inexperienced Man 


Estate planning is a job most people 
do only once. Hence they run the risk 
of making some pretty serious mis- 
takes just through lack of experience. 
Not so, however, when you have the 
helpful guidance of Fidelity’s Estate 
Planning Department. Make a date 
now—with your lawyer, your life un- 
derwriter, and Fidelity. 


FIDELITY 


PHILADELPHIA TRUST COMPANY 


Broad & Walnut Streets Philadelphia 9, Pa. 
Member Federal Deposit Insurance Corporation 


xvi 


A 
AY 
if | 
A 
3 


Journal 
The Franklin Institute 


Devoted to Science and the Mechanic Arts 


MARCH, 1959 No. 3 


Vol. 267 


THE IONOSPHERIC FARADAY EFFECT AND ITS APPLICATIONS 


BY 
FRED B. DANIELS! AND SIEGFRIED J. BAUER! 


ABSTRACT 


An outline is given of the theory of the ionospheric Faraday effect (the rotation of 
the plane of polarization of radio waves by the ionosphere). Measurements of 
ionospheric characteristics by means of lunar radio reflections and radio transmissions 
from artificial satellites utilizing this effect are discussed. The rate of fading of 
satellite signals due to the Faraday effect is found to depend upon both the integrated 
electron density up to the height of the satellite and the local electron density at the 
satellite. Thus, a possibility exists of determining both quantities provided the 
satellite transmissions are received simultaneously by two or more stations. Ex- 
pressions for the fading rate are derived and the effects of orbital parameters on the 
determination of ionospheric characteristics are discussed. 


INTRODUCTION 


Within the last few years radio reflections from the moon and, more 
recently, radio transmissions from artificial satellites have been used to 
measure ionospheric characteristics above the F-layer maximum, a 
region inaccessible to observation by vertical-incidence soundings. 
Since the Faraday effect (the rotation of the plane of polarization by 
the ionosphere) is usually made the basis of these measurements, this 
paper has been prepared for the purpose of analyzing this effect in 
detail. 


THEORY 


Magneto-Ionic Splitting 


A linearly polarized radio wave traveling through the ionosphere is 
equivalent to two components (the ordinary and extraordinary waves), 
which, in the general case, are elliptically polarized with opposite senses 


1U. S. Army Signal Research and Development Laboratory, Fort Monmouth, N. J. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
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of rotation, and which travel with different velocities. In examining 
the behavior of the two components, it is useful to distinguish between 
two cases: the longitudinal case, where the direction of propagation is 
parallel to the Earth’s magnetic field, and the transverse case, where the 
direction of propagation is perpendicular to the field. 

In the longitudinal case, if the wave frequency is sufficiently greater 
than the collision frequency of the ionospheric electrons, the major axes 
of the two ellipses are mutually perpendicular. Also, if the wave fre- 
quency is much greater than the gyromagnetic and plasma frequencies, 
the ellipses become circles, and the resultant of the two circularly 
polarized waves is linearly polarized. However, because the two com- 
ponents have different phase velocities, the plane of polarization of the 
resultant gradually rotates as the wave progresses through the iono- 
sphere. This phenomenon is known as the ionospheric Faraday effect. 

In the transverse case the two elliptical components degenerate into 
linearly polarized waves, one component always remaining parallel to 
the magnetic field and the other one perpendicular thereto. However, 
because the components travel with different velocities, the resultant 
gradually changes along the path from a linearly polarized wave to an 
elliptically polarized one, and finally to one which is circularly polarized. 
Then the sequence is reversed and the wave changes back to a linearly 
polarized one. Now, however, the plane of polarization of the resultant 
is perpendicular to the original direction. Each time that the resultant 
passes through the cycle described, there is a 90°-rotation of the linearly 
polarized resultant. 

For frequencies above 20 mc and values of less than 80° for the angle 
between the propagation direction and the magnetic field, the polariza- 
tion is nearly circular, and the Faraday rotation may be computed as 
follows: If \, and X, are the wavelengths corresponding to the two phase 
velocities v, and »v,, and dr is an element of path length, the rotations of 
the two components are given by 
oo = and = 


Figure 1 shows the relation between ¢., ¢., and d@, the angle through 


REFERENCE PLANE 


Fic. 1. Relation between phase vectors of magneto-ionic components. 


: 


March, 1959.] IoNOSPHERIC FARADAY EFFECT 189 


which the plane of polarization rotates as the wave travels the distance 
dr. From this figure, 


+ dd = — do 


do = de). 


Substituting for ¢, and ¢, and expressing \ in terms of index of 
refraction n, angular frequency w, and free space velocity c, and using 
the relations \ = v/f and n = c/v, we have: 


do = = (n, — n,)dr. (1) 


2c 


For frequencies sufficiently above the collision frequency and the 


plasma frequency (4 > 1 ) the Appleton-Hartree expression for the 
N 


indices of refraction reduces to (1)? 


No =1- i 
. 2 3 2 2 1/2 


WN 2w yn? Wy? 


where w = 2z2f, the angular frequency of the transmitted wave, 
wy = angular plasma frequency, w, = angular gyromagnetic frequency 
and 6 = angle between the direction of propagation and the Earth’s 
magnetic field. 

All practical cases can be treated by approximations of the Appleton- 
Hartree formula for the two limiting cases of longitudinal and transverse 
propagation (that is, @ = 0 or @ = 2/2). The conditions under which 
quasi-longitudinal (QZ) or quasi-transverse (Q7) propagation prevails 
are given by: 


2 
(QL) 
2 w 
and 
2 
sin tan >> (QT) 


Faraday Effect for Quast-Longitudinal Propagation 
The approximation of the Appleton-Hartree formula that holds for 
the QL case is given by: 


2 The boldface numbers in parentheses refer to the references appended to this paper. 
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= electron density (cm~) 
electronic charge (esu) 
= electronic mass (gm), and 


Earth’s magnetic field (gauss). 


Substituting in Eq. 2, we get: 


cos 6 
w* 
and for Eq. 1: 
le e*BN cos 6 


Integration of Eq. 3 results in: 


03 R 
has ge f BN cos 6dr (4) 


where gx is the total rotation along the propagation path from 0 to R. 

The above expression gives the correct result for d¢ only when the 
two components are circularly polarized. When @ is quite close to the 
limiting value that separates the QL and QT cases, the circles become 
ellipses, and when the eccentricity is large, d@ deviates from the value 
computed from Eq. 3. (In this case, d@ depends not only upon 6, but 
also upon the orientation of the plane of polarization of the incident 
wave relative to the plane containing the Earth’s magnetic field and the 
direction of propagation. ) 


Faraday Effect for Quasi-Transverse Propagation 


The approximation of the Appleton-Hartree formula that holds for 
the Q7 case is given by 
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where 


For this case, therefore, 


2 2 
WN'WH 
sin’ 6. 


No Ne 


Substituting for wy and wy gives for Eq. 1 


doar = sin? 6dr ©) 


and integration of this expression results in 


4 


e 
der = sin? 6dr (6) 


where ¢@r is the total rotation along the propagation path from 0 to R. 

The above expression for the Faraday rotation is only formally true, 
since it assumes approximately circular polarization. The actual 
situation in the case of transverse propagation has been described above 
under the heading of ‘‘Magneto-Ionic Splitting,’’ where it was pointed 
out that the polarization ellipses degenerate into lines, but the resultant 
gradually changes along the path from a line into a circle. This change 
occurs without any accompanying change in the orientation of the major 
axes of the ellipse. Therefore during this interval there is no directly 
observable change in the quantity ¢. We may, however, interpret the 
transition from linear to circular polarization as being equivalent toa 
45°-change in ¢. This interpretation is reasonable because Eq. 5 gives 


dd = 4 
from a line to a circle. Furthermore, for a dr-value twice as great, the 
circle changes back to a line that is rotated 90° with respect to the 
original line. 

The conditions under which the QL or QT approximations can be 
used have already been pointed out. It is also of interest to determine 
the relative magnitudes of the Faraday rotations for the two limiting 
cases: pure longitudinal and pure transverse propagation. 

From Eqs. 3 and 5, setting 6 = 0 and @ = 7/2, respectively, the 
ratio of dé, to ddr becomes 


for a value of dr which can be shown to correspond to the change 


do, _ 4xmcf 
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For a frequency of 20 mc, the lowest frequency transmitted by the 
Russian satellites, this ratio is approximately equal to 30. The sig- 
nificance of this result is that as the satellite approaches the position 
for which the transverse case applies, the rate of fading will drop 
sharply. This drop can be seen in the fading curves of Australian 
satellite observations (2). 


APPLICATIONS OF THE [ONOSPHERIC FARADAY EFFECT 
Moon Reflections 


Murray and Hargreaves (3) first pointed out that the slow fading of 
lunar radio echoes is caused by the Faraday effect, and that the effect 
could be utilized to measure the integrated electron content in a vertical 
column of the ionosphere. Such measurements have been made by 
Brown et al. (4) and Evans (5) on two closely spaced frequencies. 
More recently the writers (6) determined the electron content from 
direct measurements of ¢ on a single frequency. 

In moon-echo experiments, the frequencies used are usually con- 
siderably higher than the plasma frequency, the gyro frequency, and 
the collision frequency. The approximations made in deriving Eqs. 4 
and 6 are therefore valid. Furthermore, at the locations where most 
of the experimental work has been done, the magnetic dip has a value 
such that the quasi-longitudinal expression for the Faraday rotation 
applies. Inasmuch as we are interested in determining the total electron 
content in a vertical column through the ionosphere, it is convenient to 
replace the path element dr by dh sec 6, where dh is a height element, and 
6 is the zenith angle of the propagation path within the ionosphere. 
If we make this substitution, put in the numerical values of the con- 
stants, and multiply by a factor of two to allow for the round-trip path, 
Eq. 4 becomes 


@ = 4.72 X 104-f-. NB cos sec (radians) 


where the limits of integration correspond to the region in which the 
ionization is large enough to make an appreciable contribution to the 
rotation. Since B cos @ sec 6 is essentially height-independent in the first 
1000 km above the Earth’s surface, this quantity can be taken out of 
the integral.* We can therefore write 


he 
= 4.72 X 10*-f-*B cos 6 sec Ndh 
Ay 


which is the equation that is customarily used in evaluating the results 
of lunar echo experiments (5, 6). 


3 Because of the substitution the expression for ¢ will be inaccurate close to the horizon. 
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Satellite Transmissions 


For the one-way transmission from a satellite to the Earth, Eq. 4 
becomes‘ 


= 2.36 X 104-*-B cos sec Nah 
At 


where the upper limit of the integral (the height of the satellite) is now 
a function of time. To investigate the polarization fading that results 
when the signal is received with a linearly polarized antenna, we 
differentiate the above expression with respect to time. This results in 


_ wpe pl? 
2.36 X 104:f-?-B ry (cos sec 6) Ndh 


+ (cos 6 sec 5)N, | (radians/sec) (7) 


where J, is the electron density at the satellite. (In evaluating the 
derivative, the assumption is made that JN is a function of h, but not 
of t.) In this expression, the effects of horizontal and vertical motion 
are separated, and the contributions of each of these components of 
motion may be estimated if information is known regarding the orbit. 

The value of dh,/dt can be obtained from the expression for the 
radius vector of an elliptical orbit: 


a(1 — é) 


= 
1+ ecosV 


where a is the semi-major axis of the orbit, 7, is the radius vector, e 
is the eccentricity of the orbit, and W the angular variable. For 
dh,/dt we therefore have 


dh, _ ar, _ all — sin ¥ 
dt dt (1+ecosW)? dt T 


where 7 is the period of revolution. 

This result shows that dh,/dt varies sinusoidally from a value of zero 
at perigee and apogee to a maximum value of 27ea/T when ¥ = 90° 
and 270°, that is, at a height h = a(1 — &) — r., where r, is the Earth’s 
radius. (It should be noted that dh,/dt changes sign in passing through 


4 Because of the approximations made in deriving Eq. 4, ¢ will be in error by about 20 
per cent at 20 Mc (the lowest Russian satellite frequency) ; for higher frequencies this error will 
become progressively smaller. 
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perigee and apogee.) One practical conclusion that can be drawn from 
these considerations is that the second term in Eq. 7 approaches zero 
for an observer located near perigee or apogee. (It will also approach 
zero as the satellite moves outside of the ionosphere, or if the orbit is 
circular.) In the general case, however, the separation of the two terms 
in Eq. 7 may be difficult. This problem will be discussed in more detail 
below. 

Since the radial velocity component u of the satellite with respect to 
the observer introduces a Doppler frequency shift which differs for the 
two magneto-ionic components, we should investigate the contribution 
of the Doppler effect to the Faraday rotation. Because of the Doppler 
shift the two components have the frequencies 


+ 
+ ni) 


where the primed indices of refraction refer to values at the satellite’s 
location. Using these frequency values in the expression for the 
Faraday effect, we obtain 


fo 


and 


fe 


dp = — 


dr —n.) + (non. — 


(n, — Ne) (1 + 


which shows that the error due to the Doppler effect is negligible, since 


or 


2u . 
— is small compared to unity. 
c 


In order to evaluate accurately the Faraday rotation we must 
eliminate the contribution of satellite spin to the observed fading. 
This is easily done if measurements are made on two transmitted fre- 
quencies which are sufficiently separated, because the Faraday rotation 
is inversely proportional to frequency, while the contribution of the 
spin is usually the same on both frequencies (depending upon the 
antenna configuration). (This discussion assumes that the propaga- 
tion path is essentially the same at both frequencies.) 

In considering the application of the Faraday effect to the determi- 
nation of ionospheric parameters, we must examine the relative im- 
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portance of the two terms in Eq. 7. We have already pointed out that 
the second term can be neglected if dh,/dt, N., or both are small. In 
the general case, however, both terms will contribute and cos 6 sec 6 and 
its time derivative must therefore be evaluated. These quantities have 
been computed for overhead passages at Huntsville, Alabama for a 
35°-orbit (U.S. satellites) and at Belmar, N. J. for a 65°-orbit (Russian 
satellites). The results are shown in Figs. 2 to 5, inclusive. From these 
figures it is obvious that the quasi-polar orbit exhibits much greater 
values for cos 6 sec 6 and its time derivative, and is therefore the orbit 
of preference for the determination of ionospheric characteristics. For 
the quasi-polar orbit, the following cases must be considered. 


BELMAR, N. J. 
65° ORBIT N-S PASSAGE 

CURVE A:2 (cong sec 3) {Nan 
CURVE B: (cos@ sec 8 
SF CURVE C: SUM OF A ANDB 


x10*radians 


| 


NEAREST APPROACH 


5 4 5 6 

TIME (MINUTES) 
Fic. 6. Computed fading rate for NW-SE satellite passage over Belmar, N. J. 
as a function of time. 


Case I-—Satellite well above the I-layer maximum at the point of 
observation. We have already pointed out that for this case, N, is 
small, and the main contribution to d¢/dt in Eq. 7 comes from the term 
which contains {Ndh. 

Case II--Satellite well below the F-layer maximum at the point of 
observation. For this case, both NV, and /Ndh are small; the total 
Faraday effect will be small, and this case is therefore of little practical 
value in determining ionospheric parameters. 

Case III—Satellite near the F-layer maximum at the point of ob- 
servation. For this case, the terms in Eq. 7 have been computed as a 
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function of time for an overhead passage at Belmar, N. J. The follow- 
ing conditions have been assumed: North to South passage, h, = 350 
km, dh,/dt = 5 X 104 cm/sec (a representative value for the orbit of 
Sputnik II), N, = 7 X 10° electrons cm-*, and f{Ndh = 5 X 10" 
electrons cm~*. The results of this computation are shown in Fig. 6 
in which the fading rate ¢ has been normalized with respect to frequency. 
Curve A shows the contribution of the term containing /Ndh, and 
Curve B the contribution of the term containing N,. It is obvious 
that the two terms, for the case assumed, are of approximately equal 
magnitude, and that the computed Faraday fading (Curve C) is the 
combined effect of these terms. It should be noted that observed 
Faraday fading rates for satellite passages above the F peak closely 
resemble Curve C (7) (see Fig. 7). 
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20Mc/s FADES PER MINUTE 


0332 334 0336 0338 0340 
GMT i9TH OCTOBER 1957 

Fic. 7. Observed fading rate due to Faraday rotation. (Reproduced, with permission 
of the Institution of Electrical Engineers, from B. Burgess, Proc. Instn. E. E., Vol. 105, Part B, 


p. 113 (1958).) 


kor Case I, the variation in fading rate will be given by a curve 
having the same general shape as Curve A of Fig. 6. As far as Case II 
is concerned it has already been pointed out that the Faraday rotation 
will be too small to be of practical use. 

One possibility exists, at least theoretically, of separating the two 
ionospheric quantities /Ndh and N,, namely, to observe the satellite 
transmissions simultaneously at two stations separated by a sufficient 
distance. This method is only applicable when the satellite is ob- 
served at a point where dh,/dt is close to its maximum value, and 
furthermore, would not be applicable for a circular orbit for which 
dh,/dt = 0. 
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The Faraday fading rate observed at stations 1 and 2 can be ex- 
pressed by 


k(ai + 5,N,) 
+ 62N,) 


$1(t) 
$2(t) 


k = 2.36 X 10*f-?-B 


0 
(cos sec 6) 


b 


(dh,/dt) cos @ sec 6. 


Solving these equations leads to 


_ — 


and 


In practice, this method will only be feasible as a means of measuring 
local electron densities well above the F-layer maximum if the dh,/dt 
of the satellite is large throughout the region of interest. 

The above discussion applies to the case where the satellite trans- 
missions are received with a fixed, linearly polarized antenna. Another 
method of utilizing the Faraday effect has been proposed by one of the 
authors (8). This method is to transmit two closely-spaced frequencies 
and to receive each signal with a pair of circularly polarized antennas 
having a continuously variable phase separation, so as to be the equiva- 
lent of a continuously rotating dipole. From the time separation of the 
signal minima on the two frequencies, the integrated electron density 
may be computed. This method has some advantages, but has the 
obvious disadvantage of requiring complicated specially designed re- 
ceiving equipment, whereas the method described in the present paper 
requires only simple dipole antennas and standard communications 
receivers. 


CONCLUSIONS 


The Faraday effect is a very useful tool for the determination of 
ionospheric parameters that are normally not obtainable by other 
methods. For the determination of the integrated electron content we 
may use either lunar radio reflections or transmissions from an artificial 
satellite. The satellite, however, would have to be placed in an orbit 
all parts of which lie well above the F region, if the total electron content 
of the ionosphere is to be determined on a world-wide basis. 
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The local electron density can be determined from satellite trans- 
missions, but the satellite orbit must extend from below the F-layer peak 
to well above it, so that dh,/dt is large throughout the region of interest. 
In addition, two or more appropriately chosen stations must receive 
the satellite signals simultaneously to allow separation of the inte- 
grated and the local electron densities. Maximum Faraday rotation 
will be obtained with polar or near-polar orbits, for propagation direc- 
tions that are not too nearly perpendicular to the Earth’s magnetic 
field, and for frequencies that are low (because of the f-? dependence) 
but sufficiently above the plasma frequency of the F layer. 
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ON EXPLANATIONS OF ELECTRIC AND MAGNETIC 
CONSTANTS AND UNITS 


BY 
A. T. GRESKY! 


ABSTRACT 


Employing dimensional methods, it is suggested that “explanations” of the free- 
space dielectric constant e and the magnetic permeability constant » may lie in their 
proposed relationships to four highly fundamental quantities; for example, e, = m2 
(hac) and py = ha,(m2c)“. The latter is considered to reflect a specific intrinsic 
property of the electron, that is, a constancy in the ratio of its magnetic moment (Bohr 
magneton) to its spin (4/2). The former, equal to 1/y,c?, reflects an additional intrin- 
sic relationship of the electron to light velocity, presumably a characteristic of its 
associated photons. Other constants and units are similarly analyzed to suggest pos- 
sible methods for understanding the submicroscopic nature of electric and magnetic 
phenomena. 


INTRODUCTION 


The mysteries of interaction-at-a-distance, as encountered in gravi- 
tation, electrostatics, magnetostatics, or nuclear binding, are major 
items of interest to scientific philosophers. The search for general laws 
of physics and unified field theories inevitably must deal with these 
mysteries, which scientists usually reconcile or rationalize by use of the 
mathematical imagination. Explanations of the interaction phenom- 
ena are considered to reside in certain unexplainable free-space proper- 
ties, as exemplified in constants such as G, e, and p. 

Whereas the constant G dealt with forces between particle masses, 
and thus acquired a knowledgeable value and dimensions of 6.67 -10~* 
dyne-cm?-g~?, the constants « and yu, dealing with forces between ab- 
stractions like charges or magnet poles, were necessarily assigned values 
of dimensionless unity. Proper numerical coefficients and dimensions 
could then be arbitrarily assigned to units of electrostatic charge (stat- 
coulomb Q = 1.0e!/*g'/?cm*/*sec-!) or magnetostatic poles (maxwell 
m = 1.0 p!/2g!/2cm3/*sec—!) and thus arrive at proper dimensions of unit 
force (dyne) from defining expressions like Q,Q2(er?)~! or mym2(ur?)-?. 


DISCUSSION 


When measurements of the electronic mass m, and the e/m, ratios 
ultimately established that the individual electron’s charge e was equal 
to 4.8-10- statcoulomb (and/or that one statcoulomb involved the 
combined charges of N., = 2.08-10° electrons), there were some in- 
vestigators who speculated that the electrostatic equation could have 
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been alternately written as 


(Nal) (Net): - (Name) - 1.02 

er er? "sec? 
Whereas conventionally « = 1.0, it would then be seen that e¢, in the last 
expression, where 7 = 1.0 cm, would necessarily acquire a coefficient 
and dimensions equal to 3.595-10-** g sectcm-*. Noting that ¢,~', 
equal to 2.77-10** cm*g~'sec~*, had dimensions identical to G = 6.67 
-10-* cm*g~'sec~2, those who had already recognized the similarity of 
the equations for electrostatic and gravitational forces, were struck with 
the possibility that the former as well as the latter dealt with mass 
interaction-at-a-distance. However, this concept was strongly dis- 
couraged by the fact that (€,(7)~', equal to 4.15-10*, denoted a fantastic 
deviation from experience in models of gravitational forces. 

As suggested in a previous paper (1)?, if specific physical laws in 
classical mechanics were granted to be entirely different in nature from 
the specific laws of electromagnetism or quantum mechanics, the vari- 
ance of 10 need not be an absolute deterrent to philosophic speculation. 
When it is recognized that the quantity 


g sec? 


= = 3.5 —36 
= = 3.5048-10-1 


of 


(signifying that e, is peculiarly related to the very fundamental quanti- 
ties of electronic mass m,, angular momentum ih, the fine structure con- 
stant a,, and light velocity c), one might reasonably ask whether the 
dielectric property of free-space cannot be perspicuously explained in 
relatively simple terms of mass-space-time ; that is, in a similar sense to 
the “explanation” of gravitational properties by the use of G. 

It was shown (1) that G, formulated as \*vo2M~—', can be considered 
as a simple, quantitative, and dimensional statement of the harmonic 
law of astronomy, indicating a mass-space-time~! relationship of a 
central massive body to the orbits of its less massive satellites. A new 
“astronomic”’ constant pw, = 1.2948-10-8 cm formulated as 
(Myv,)~-!, was also postulated (1) as a similar modified statement of 
Blackett’s hypothesis (4,5) which related a solid body’s rotational 
motion to its polar magnetic field intensity. Its dimensions of cm g~', 
or cm sec(g sec)~!, were indications of yet another quantitative mass- 
space-time! law of macroscopic bodies in free-space. Since the mag- 
netic field intensity v,, could be expressed in sec~! dimensions (gauss 
= 1.759-107 sec!) it was possible to visualize and interpret a new di- 
mension of time (or time~!) in the measurement of massive particle 
properties and mechanics. 


? Boldface numbers in parentheses refer to the references appended to this paper. 
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That the gauss unit could be expressed in sec~! dimensions resulted 
from a study of ‘‘mechanical equivalents” for units and constants em- 
ployed in electricity, electromagnetism, electrostatics, and magneto- 
statics. This study revealed, as for ¢, above, that many of these equiva- 
lents could be related to m., a., #, c, and numbers (N.., Na, Nem) known 
from em,.-! measurements, for example : 


specific magnet pole 8.44915-10-8 
Me sec 


2 
sec 


maxwell (unit pole) 


e 


cm 
(of vacuum) 320955: 10" 


e 


maxwell-cm Nehacm _ 594.107 


e 
ha, cm? 


Bohr magneton = 1632-10-13 — 


mec sec 


Nuha, 1.7594-107 


auss 
8 mcm? sec 


abamp Name _ 5.6837-10-*-£., 
cm sec 


oersted 
cm? cm sec 


statcoulomb N.m, = 1.896-10-'8 g, 
abcoulomb Nenm, = 5.6837-10-8 g, 
coulomb Nam, = 5.6837-10- g, ete. 


The observed relationship between uw, and the Bohr magneton sug- 
gested that the former was a measure of the ratio of an electron’s mag- 
netic moment to its spin .S, (or #/2) and that the 


“m3 
Bohr magneton = y,S, = = 1,632-10- 
mM sec 


Since several investigators (2-5) have considered that the magnetic 
field intensities of both macroscopic and submicroscopic particles are 
proportional to their mass and angular momentum, it was finally con- 
sidered that 
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sometimes (6) called the electron’s classical radius, might actually be a 
universal constant of submicroscopic particle physics relating mass 
‘magnetic moment to particle spin. This indicated the possibility that 
yet another quantity M, = m.a,.“', evident in L,,, was also a universal 
constant relating any submicroscopic particle’s mass to its fine structure 
constant. It was thus conceived that a, is a specific electronic quantity 
which might define a constancy in the relationship of the electron’s 
magnetic field intensity v, to its angular velocity », (for example, as 
surmised from », = E,h-"), that is, 


Me Van 
= = — = e —3 
= ( ), 7.2978-10-%. 


Another quantity, Fy = Gu,-* = 3.9782-108 g cm sec~’, relating pw, to 
the gravitational constant G, was considered (1) to be a general constant 
of physics, since, also, 


suggesting the existence of a single ‘“‘law’’ of quantum mechanics, electro- 
magnetism, and classical mechanics. (It has been noted (7) in the press 
that Dr. Werner Heisenberg’s unified field equation might depend im- 
portantly on three constants: that is, #, c, and an unknown constant x, 
with dimensions of length approximating the 10-"* cm range of nuclear 
forces. Since the above quantity L,, = 2.8185-10-' cm has properties 
suggestive of the constant x, and since # and c also appear in quantity 
F,, one may speculate whether the latter is an essential or component 
part of Dr. Heisenberg’s hypothesis of the unified field.) Returning to 
the “‘explanation”’ of the dielectric constant e, and the magnetic per- 
meability constant yy, it is of interest to note also that 


3.5948 - 10- 


Lm _ 3.0955-10"°™, and 
m, g 


_ here, 


m,? 


With knowledge of L,,’s meaning, it may then be surmized that e, re- 
flects a fundamental relationship between the electron’s magnetic mo- 
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ment, spin, and light velocity (presumably associated with a character- 
istic photon of rest mass m, with energy m,c? equal m,v*?). Similarly, 
u, would appear to reflect a fundamental relationship between the elec- 
tron’s magnetic moment and its spin. The general constant F, may 
then be shown to relate an electron’s energy E,(m,v? or m,c*), a wave- 
“radius” Xo (characteristic of its energy, according to acE,-'), its mass 
m,, its magnetic moment M,, and its spin S,; that is, 


_ EXoS2 _ energy: wave-radius- spin? 


m?M,? mass?- magnetic moment? sec?” 


The several interrelationships observed among (1) the general con- 
stants #, c, M,, Lm, F,, (2) the electronic constants m,, a., S., Ms, ug, 
and e,, and (3) the numbers N,,, Nem, and N,., may possibly afford clues 
to “explanations” of all electric and magnetic phenomena. A sug- 
gested approach to the use of such clues is outlined in the Appendix. 
Although the highly developed technology of electronics and mag- 
netism has not required actual explanations of space properties such as 
€, w, Or resistivity, it may be that electromagnetic science and theory 
would profit from any possible resolution of their mysteries. 
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APPENDIX 
Hypothetical Models of Electron Mechanics and Magnetism 


Quantities (1) through (8) in Table I may be suggested to describe the rotating electron 
at three different energies and/or voltages (abvolt, volt, and 300 volts); that is, characteristic 
of potential differences or electromotive forces encountered in the unit systems of electromag- 


TABLE I.—Mass, Space, and Time™ Properties of Electrons at 10-8, 1.0, and 300 Ev 


10-8Ev Electron 1.0Ev* Electron 300Ev Electron 
(Emu System) (Pract. Elect.) (Esu System) 


(1) m, = electron rest mass, = 9.11-107%8 9.11-10-%8 10-*8 g 

(2) Ae = A(2m,c) = S,(mc) = 1.93-10-" 1.93-107! cm 

(3) Ar = (h?/m,E,)'? = 2.76-10~4 2.76-10-8 cm 

(4) » = Eh = 1.52. 1.52-10'5 -10!7 sec! 

(5) om = = 1.11- 1.11-10'8 

(6) = Eh = Eh 1.52- 1,52-10'5 -10!7 sec! 

(7) Xo = AcE," = 1.97- 1.97-10-5 -10-§ cm 

(8) m, = E,c? = 1.78- 1.78 - 107% -10-% g 

(9) = E, = 1.60- 1.60-10712* -10- g cm*sec™* 
(10) = E, = 1.60- 1.60-1072* -10-° g cm*sec™ 
(11) Aw = 3.00- 3.00 - 10" .00- 10! cm sec™ 
(12) = = h = 1.05- 1.05- 10727 .05- 10-27 g cm?sec™! 
(13) md.c = 2/2 = S, = 5.27: 5.27- .27-10-*8 g cm*sec™! 
(14) AZ», = = 1.16- 1.16-10° -10° cm?*sec™ 
(15) = hams! = Ln 8.45-10° 8.45-1078 cm’sec™ 
(16) = = Lin 2.82. 2.82 -10-8 82-10-83 cm 

(17) APAa = = 2.94- 2.94- 10-26 .77 +107 cm? 

(18) = My 1.63- 1.63-10-% cm'sec™ 
(19) = a, = 7.29. 7.29-10-3 .29-10-% 

(20) md’vm?/2 = Enya. = 2.18- 2.18-10-"! -10-" cm*sec™ 
(21) = wg = 3.09 - 10" -10'* cm go 
(22) = € = -10- 3.59 - 10736 -10-%6 g sectcm™? 
(23) N = no. of electrons. 6.24-10'8 -10° 

(24) NE, = NE, = E, = 1.00-107 .00- 10° g cm*sec™* 
(25) A277 = Em, = V, = 1.76-10'5 10'7 cm*sec™? 


* One electron-volt, equal to 1.6-107"? erg. 


netism, practical electricity, and electrostatics. The rest mass m, is considered to be located 
in a dense spherical “nucleus” of \q radius, shown as constant for all energies E,. Note that 
4rd, would be equal to the Compton wave-length \ = h/m,c, and that dq is equal to #/2 (half- 
integral “spin’’) divided by mc, suggesting a relationship of “spin’’ to m,’s original creation 
near light velocity. The “equatorial” or effective radius \,, shown to increase with decreasing 
energy, will be noted as equal to the de Broglie wave-length, \ = h/my, divided by 27, Quan- 
tity (17) suggests that m,’s effective volume, dependent on the product of \g (axial radius) and 
4,2, becomes increasingly oblate with decreasing energy. The value for », (angular velocity, 
in radians-sec™) will be observed as 27 times greater than the wave-frequency », equal to E,h™, 
which would characterize the given electron as a wave of radiation. Quantity (8) suggests that 
characteristic photons of rest mass m,, with energy E, equal E,, become associated with elec- 
trons at all energies. (Note that m, + m, does not become appreciably greater than m, until 
the energy approaches m,c?, at which point m, + m, would approach 2m,.) The photon of rest 
” or standing-wave with angular velocity vo 
equal to m,’s angular velocity », and a radius Xo equal to cvo or AcE. (It may also be 
imagined as revolving about m, in an orbit of radius A» with angular velocity vo, exhibiting a 


mass m, may be imagined as a rotating “particle 
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peripheral velocity \ovo equal toc. If emitted as radiation, it would appear as a light particle- 
wave of frequency vo(27)~ and length 27Ao.) Finally, m,. may be considered to possess a mag- 
netic field, polarized along axis \g and perpendicular to the equatorial plane A,*, having an 
intensity vm equal to v,a, vibrations (?) per second. The measure of this field, in gauss, for the 
three given electrons may be found by dividing the v, values by 1.759-107 sec™! gauss”. 
Note, from quantities (15) and (18), that each free electron may then be characterized as a 
magnetic pole equal to 4.8-10~- maxwell and as having magnetic moment of 9.3- 107?! max- 
well-cm, or one Bohr magneton (where, as previously (1) suggested, the maxwell and the Bohr 
magneton are equivalent to 1.759- 107 cm*sec™! and cmsec™, respectively). 

By formulation, quantities (12) and (13) suggest the similarities that exist in the constants 
of angular momentum & and spin S,. Quantities (19) through (22) suggest the formulation of 
electronic constants such as the fine structure constant, the Rydberg energy, and the magnetic 
permeability and dielectric constants of free-space. Note that quantity (16), sometimes (6) 
referred to as the ‘‘classical electron radius,"’ may appear to be a universal constant relating the 
mass, magnetic moments, and spins of electrons (and perhaps of other submicroscopic particles). 

Quantity (23) may be determined from knowledge of charge-to-mass ratios, and quantity 
(24) defines the energy units (erg and joule) in the cgs systems of measurement. Quantity 
(25) suggests the numerical and dimensional values of mechanical equivalents for the familiar 
abvolt, volt, and statvolt of potential difference or electromotive force. 


Suggested Definitions of Units and Constants of Electricity and Magnetism 


An understanding of electron properties as shown in (1) through (8) above, as well as cer- 
tain interrelationships as suggested in (9) through (25), may be helpful in interpreting the 
fundamental implications of the suggested (1) mechanical equivalents. For example, where 


me = 9.1055-10~*8 g, 


Nes = 2.0823-10°, 


Nem = 6.242-101, a, = 7.2978-10-3, 

No = 1.0, h = 1,0542-10-?7 g cm’sec™!, 

r = 1.0 cm, c = 2.9977-10' cm 

= 1.0 cm, M, = ma," = 1.2477-10-* g 

a =1.0cm, Lin = ha,(m,c)™ = 2.8185-10-8 cm, 
v = rt = 1.0 cm sec", t = 1.0 sec, 


and where r, a, and / represent mutually perpendicular directions of separation or length, one 
may alternately derive equivalents for quantities employed in electrostatics, magnetostatics, 
and electromagnetism, as shown below. 


Electrostatic Quantities 
+ + 2 
Qes Qes (Neme ) ) het Ne hei (26) 


r cm? “Nom? Nocm? 


suggests an equivalent of 1.896-10~"8 g for the statcoulomb, 1.0 g cm sec? for the dyne of elec- 
trostatic force, and 3.5948 - g for the dielectric constant of free-space. Note, in 
this case, that Q.,~ refers to a number of free electron masses, while Q..* may refer to a similar 
number of positrons (of rest mass m,) located near the nuclei of ionized atoms. Note that ¢, 
might also have been written as Nom.M,(hc)™ or Nom,(L»c?), and that dielectric values 
greater than ¢, (vac.) would suggest an increase in the number No. Note that an equivalent 
unit of potential energy or work would be equal to F,,r or one erg, and a unit of surface charge 
density might be written as Q,,(al)™, or Nesmecm™?, or 1.896-10°'8 g The unit of electric 
field intensity A, might be written, as follows: 


N.am,* 
cm? 


x Ime? _ Ne 


2 
Nom.  Nocm? Lmc?, (27) 


bi 

eS 

1 

Qua 4 
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suggesting an equivalent of 5.2742-10!7 cm sec~ (dimensions of acceleration) for the statvolt- 
cm~, The unit of potential difference V,, the statvolt, would be equal to Agr or 5.2742-10"" 
cm*sec™, and the unit of capacitance C, would be equal to 


1 Nocm Nom cm 


suggesting an equivalent of 3.5948-10-*6 g sec*cm™ for the statfarad unit. 
Magnetostatic Quantities 


a? Mg cm? NoLm Nocm? 


m,-m,* x 1 os (NeLme)~ Linc? (29) 


suggests an equivalent of 1.7594-107 cm*sec™ for the unit magnetic pole or maxwell, 3.0955 - 10" 
cm g™ for the magnetic permeability constant yw, of free-space, and 1.0 g cm sec~ for the dyne 
of magnetostatic force. Note that the distance of pole separation is symbolized by a rather than 
r as in the electrostatic case, and that the unit of magnetic field intensity B, would be written 
as 
mo _ _ (30) 
rl rl cm? 
suggesting an equivalent of 1.7597-107 sec™! for the gauss, or maxwell-per-unit-area. Note 
that uz, might also be written as Noha.(mc)™, and that it is equal to 1/ec*. Also note, for 
permeabilities greater than yw, (vac.), that an increase in the number No would be indicated. 
The maxwell-cm of magnetic moment MM, could be written as 


ma = B,arl = M, = NuLmca (31) 


suggesting an equivalent of 1.7594-107 cm*sec™, equal to 1.08-10” Bohr magnetons. 


Electromagnetic Quantities 


cm? _ Neale 


Baul rl Voi cm? sec sec 


suggests an equivalent of 1.7594-107 cm*sec™? for the abvolt unit of electromotive force. Note 
the symbolic indication that such a force is produced along the length / of a wire moving with 
velocity rt (in the r direction and perpendicular to magnetic lines of force oriented in an a 
direction). The equivalent for the abvolt-cm™ unit of electric field intensity A, along / would 
be or 1.7594-107 cm sec™. 


Ln es 
Nolom Neate (32) 


B, = Bolly; Nor? 


suggests an equivalent of 5.6837-10-* g cm™ sec! for the oersted unit of magnetizing force H,, 
where a unit current 
Hor, Name _ Name ., 


equivalent to 5.6837 -10-8 g sec"! or one abamp, flowing in a conductor arc of length /, produces 


a field of gauss intensity at radius r, ina medium of permeability uz. Note that J, can also be 
written as Nemm,sec™. From 


_ Neamt y, None (34) 


sec Nem™, sec 


V, = 


an equivalent for the abohm of resistance R, would be equal to 3.0955-10' cm? g“ sec“. Note 
that the abohm-cm of resistivity R¢ might then be written as NesLmcl/Nemm. (ar) or 3.0955 
-10 cm*gsec™, implying that it measures an empirical ratio of electron magnet poles per 
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unit length to electrons per unit area of the conductor. Electromagnetic force F., between two 
parallel current-carrying wires, in vacuo, can be expressed as 


(Io!) Nedmac? Nohew _ 


mec Nor? * 


(35) 


The dyne unit of electrical force F,; can be expressed as 


we sec cm sec 


Interpretations from Mechanical Equivalences 


From examples in the above discussions, and the data in Table I, it is possible to evolve 
certain imaginative interpretations of phenomena. One example may be cogsidered, as fol- 
lows. From quantities (2), (3), (15), (29), (30), and (31), one might surtifise that: (a) a 
magnetic pole face of area rl, with gauss intensity, is characterized by 2.0823-10® polarized 
electrons (and/or electron magnet poles \,?vm) per unit area, with their equatorial \,? planes 
parallel to the face; (b) since each electron would thus occupy a maximum of 4.8-107~' cm*, its, 
minimum allowed energy might appear to be in the order of 10~7 g cm*sec™? or 6.25-10~6 ev 
that is, where its A, would not exceed about 10~§ cm; (c) if the magnet were a cube with 1.0-cm 
sides with a magnetic moment of one maxwell-cm, the foregoing would imply that cm-A,g™ or 
5.18-10" electron poles along the magnet’s length a would be required to account for 1.08- 10” 
Bohr magnetons within the cube. Thus, one might visualize a magnetic field of gauss intensity 
as containing 2.082-10° parallel “‘lines’’ of electrons per cm? of rl area. The individual lines, 
polarized along direction a and perpendicular to 7] would then appear to consist of some 5.18 
- 10 axially-connected electrons per cm of line-length. A 1000-gauss field might thus be char- 
acterized by 2.082-10" such “‘lines” per cm? of field area. 

Similarly, from (34) and its discussion, a cubic conductor with 1.0-cm sides, having a 
hypothetical resistivity of 1.0 abohm-cm (probably characteristic of supra-conductors), might 
be visualized to contain 6.242-10" electrons per cm? of face area ar, oriented with their \,* 
planes perpendicular to the face; 2.082-10° electron lines, with their axes perpendicular to the 
direction of current flow, would then be located along length /. (Note that good metal con- 
ductors have resistivities in the order of 1.0 microhm-cm, or 10? abohm-cm, indicating a less 
favorable, or higher, ratio than NesNem™.) 

From (35), concerning the fields and forces of two current-carrying wires in vacuo, an 
example of oppositely directed currents provides a mechanical concept of geometrical and polar 
relations among the three fields involved in simultaneous electrical, electromagnetic, and di- 
electric forces. The two electric fields in opposite directions / contribute to the induction of a 
magnetic field in the space between them, oriented in direction a. The magnetic fields circling 
each wire can be visualized as lines of axially-connected electrons, with line-density increasing 
near the wires, and being enhanced in the region between the wires. Relative polarization in 
the latter region would be such as to suggest magnetic dipolar repulsion of the lines, and thus a 
mechanical force on the wires. A dielectric field (perpendicular to both the electric and mag- - 
netic field directions) also exists in the r direction between the wires, however its energy would 
depend on V,?C, and would become significant only at higher voltages. 


Relating Constants and Units to the Universal Constant F, 


Since each of the various constants mentioned in the above discussions is related to m,, 
c, a, and/or h, it will be noted that each can be related to the other by use of a general equation 
exemplified below: 


3.978210 
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It will also be recognized that certain interrelationships of the various units and constants 
will also be equal to F,, for example: 


_ (abohm) (abamp) (Ao) em (oersted) = 3.9782-108 
MeptgB (electronic mass,) (ug) (gauss) i sec? 


F, 


— _ (abohm) (abwatt) (Ao)em(statcoulomb)? _ 3,9782-108 gem 


m,m,* em (maxwell)? 


LoRaFemdotg _ (abamp) (abohm -cm) (dyne) (Ao)em _ 3.9782-108 
QemMaptg (abcoulomb) (m,)em (ug) sec? 


F, 


It will be noted that the latter expression contains a relationship among the three elusive 
properties of space (or media), that is, resistivity (Rg), dielectric (€,), and magnetic permeabil- 
ity (u,). Speculatively, it may be hoped that such knowledgable relationships can eventually 
lead to a clearer perception into the nature of medial properties and the three related forces 
(electric, electromagnetic and electrostatic) which characterize electron motions and inter- 


actions. 


THE INFLUENCE OF GRAVITY ON THE FOLDING OF 
A LAYERED VISCOELASTIC MEDIUM 
UNDER COMPRESSION 


BY 
M. A. BIOT! 


ABSTRACT 


The theory of folding of a layered viscoelastic medium under a horizontal com- 
pression parallel with the layer, as developed in previous publications (1, 2, 3),? is ex- 
tended to take into account the influence of gravity. Characteristic features of the 
folding are evaluated in terms of nondimensional parameters. The influence of 
gravity appears through a parameter which depends on the relative magnitude of the 
gravity forces and the horizontal compression. Specific cases discussed are that of a 
layer on top of a half-space and the layer embedded between two half-spaces, with 
the bottom medium either denser or less dense than the top one. 


1, INTRODUCTION 


In references 1 and 2 we have investigated the folding due to in- 
stability when a layered viscoelastic medium is subject to a compression 
in a direction parallel with the layer. The layer was assumed to be 
embedded in an indefinite solid or lying on the surface of a half-space. 
The problem was analyzed from a very general viewpoint and the visco- 
elastic and hereditary characteristics of the material were assumed to 
obey the most general linear laws. In particular, conclusions are drawn 
for materials which obey the Onsager relations and the laws of the 
thermodynamics of irreversible processes as formulated in some our 
previous work (7, 8). The folding was also evaluated for various 
specific combinations of viscoelastic media, and attention was given to 
the effect of interfacial friction at the boundary and to the evaluation of 
the magnitude of the folding as distinct from mere instability (2). In 
the cited work we used an approximation of the ‘‘plate theory” type to 
represent the folding of the layer. In reference 3 this approximation 
was discarded and the problem was treated by using exact equations of 
the stability of a prestressed continuum as developed in earlier papers 
(4, 5, 6, 7). The problem of stability of the surface of a viscoelastic 
half-space was also treated as a particular case. 

In the present analysis we are dealing with the problem of folding 
of a layer, either lying on top of a viscoelastic half-space, or embedded 
between two such half-spaces. When in addition to a compression 
parallel with the direction of the layer there is also a uniform gravity 
field perpendicular to this direction, we are therefore dealing with in- 


1 Consultant, Shell Development Company, Houston, Texas. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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stability under a combined state of prestress—the horizontal compres- 
sion in the layer, and a prestress due to gravity which is assumed to be a 
hydrostatic field. The addition of gravity to the problem is essential 
of course if we have in mind some of the geological applications. We 
shall discuss these at a later date. 

A general formulation of the problem is outlined in Section 2. 
From the previously established theory (4, 5, 6, 7), we derive the equa- 
tions for the stability problem of a continuum under a uniform compres- 
sion and a hydrostatic stress gradient due to gravity. These equations 
are quite general and may be applied to both compressible and incom- 
pressible media. For simplicity we have applied the theory only to 
incompressible media and in Section 3 we show that in this case simpli- 
fications may be introduced along with a “plate type” equation for the 
layer. ‘This leads to a very simple equation for the folding of a layer 
lying on top of a half-space. The case of a layer and underlying medium 
which are incompressible and purely viscous is discussed in Section 4. 
The dominant wave length (that is, that whose amplitude grows at the 
fastest rate) is evaluated as a function of the ratio of the viscosity coeffi- 
cients of the two media and a nondimensional parameter which repre- 
sents the influence of gravity. This parameter is a ratio of two quanti- 
ties, one is the pressure of a column of underlying material of height 
equal to the layer thickness and the other the compressive horizontal 
force. The density of the layer itself does not enter into the picture 
because of our assumption that the layer behaves like a thin plate of 
constant thickness. The relative amplitude of the folding for the 
dominant wave length is also evaluated. This yields an indication of 
the magnitude of the instability. 

Section 5 discusses the case of a layer embedded between two differ- 
ent materials. The materials are again assumed incompressible and 
purely viscous. Two cases must be considered. In the first case the 
underlying material is denser and the effect of gravity is stabilizing. 
The theory is the same as for the layer lying on top of the half-space 
except for modified parameters which introduce the sum of the vis- 
cosities of top and bottom media and the difference of their densities. 
The effect of gravity tends to shorten the dominant wave length. The 
second case is that of an inverse density gradient, that is, the under- 
lying material is lighter than the top one. The layer is then unstable 
because of gravity alone. The phenomenon then becomes somewhat 
more complex. For certain ranges of the parameters the combined 
effect of gravity and compression will increase the dominant wave 
length while in another range the dominant wave length will be con- 
trolled by geometric factors such as the depth of the surrounding me- 
dium. It should be kept in mind that we have not taken into considera- 
tion any inertia forces. The deformations are therefore assumed to be 
sufficiently slow for these forces to be negligible. 


. 
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Strictly speaking, the theory is also restricted to infinitesimal 
deformations but in practice it is of course applicable to a large category 
of observable phenomena in the same sense as this is understood to be 
true in the applications of the classical linear Theory of Elasticity. 
Furthermore it is most likely that the folding wave length which ap- 
pears in the incipient phase of the instability will also correspond to the 
characteristic feature in the later phases with large deformations. An 
experimental program of model tests on the folding of layered visco- 
elastic media has been initiated. Results are in good agreement with 
the theory. They will be presented in later publications. 


2. GENERAL EQUATIONS FOR THE DEFORMATION OF A PRESTRESSED MEDIUM 


In reference 3 we have analyzed the folding instability under com- 
pression of a layered viscoelastic medium by using the exact equations 
for the deformation of an anelastic continuum under prestress. These 
equations are the same as those which we established in references 4, 5 
and 6 for the elastic continuum. The general applicability of these 
equations to anelastic media was discussed in reference 7. A brief re- 
derivation of the equations was also presented in reference 3. In the 
previous work we did not introduce the gravity force. In the following 
we shall consider the prestress to result from both the lateral compres- 
sive forces on the layered material and the hydrostatic pressure gradient 
due to its weight. In the present section we shall write out the equa- 
tions and boundary conditions for this particular state of prestress. 
We shall also restrict ourselves to two-dimensional deformations. As 
established in the work quoted above the incremental stress field of the 
continuum satisfies the equations 


OS 22 dw dw 
+ + pYw — 2S 12 ax + (Sir S22) ay 


and another similar equation obtained by permuting x and y, 1 and 2 
and changing w into — w. The coordinate system is assumed to be 
clockwise. In these equations S1;, S22, Si. are the initial stress com- 
ponents in the x, y plane, 511, S22, Siz are the incremental stresses, and 
pY is the y component of the body force per unit volume (p mass den- 
sity). The displacement field u, v yields the strain components 


Ou 
ax Cy, = (2) 


1fov a 
= 5 


2\dx dy 


’ 

om 

i 

| 

| 

A 

) 
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1 f ov Ou 
3(2-). (3) 


In this representation the stress components s;; are the incremental 
stresses at the displaced point x + u, y + v, and relative to axes rotated 
by an angle w. The strain components e;; are also referred to the same 
displaced and rotated axes. The latter are of course not distinct, in the 
first order, from the classical components for small strain. 

Consider now a particular case of prestress, namely a combination 
of a horizontal compressive stress P and a hydrostatic stress due to 
gravity. The gravity force is assumed parallel with the y axis. In that 


case 
(4) 


and the rotation 


Y=g 


where g is acceleration of gravity and the initial state of stress is 


Si = — P — pgy 
Soe = — pgy (S) 
Sie = 0. 


The origin of the y direction is left arbitrary which is equivalent to 
adding an arbitrary uniform hydrostatic stress in order to suit the 
boundary conditions. Equations 1 to be satisfied by the incremental 
stress field then become 


OSi1 Ov Ow 
— P—=0 
Ox + oy + Ox oy 6) 
OS OSo9 Ou dw 
Ox oy Ox Ox 


The boundary conditions may be formulated using the following rela- 
tions (4, 5, 6): 
dF, = — [Sie + Siz — Seow + — jdx 

+- [sir + Si; — Siew + — 
[See + Soo + Siew + — 

+ [Sie + Sie + Siw + Sirlyy = 


(7) 


dF, 


where dF, and dF, are the forces acting on an element of boundary 
initially represented by dx, dy. The force dF is acting on matter lying 
on the left side of the arc dx, dy, considered as a vector differential. 
For the particular prestress field under consideration, Eq. 5, these rela- 
tions become 
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dF, + Pes + 2° | ax 


+ E — (P + pgy) (1 + mele (8) 


dF, = — | — (1 +) + [su - Po + 5, 


We must also express relations between the incremental stresses and 
deformations. If we assume that they are of the same form as for an 
unstressed medium, we write 


= 20e:;;+ 6:;Re (e = + Cyy) (9) 


with the operators (p = 5 , time deriv ative), 


R= + R + R’p. 


We have previously derived these expressions from irreversible thermo- 
dynamics (8). Strictly speaking their use for prestressed media may 
involve an approximation. ‘The nature of the limitations involved was 
discussed in reference 2... However, in most problems these limitations 
are only of academic interest. 

Substitution of expression 9 for the stresses s;; into Eqs. 6 yields two 
operational equations for the displacements u and v. The problem of 
deformation and folding of a layer lying on a semi-infinite medium both 
of heavy viscoelastic material and subject to a horizontal compression 
is thus reduced to solving the two simultaneous equations for u and v 
with appropriate boundary conditions. This can be done in a way 
entirely analogous to the procedure followed in reference 3. The exact 
solution however is somewhat involved and we have preferred to follow 
an approximate but much simpler approach in analogy with our treat- 
ment of the problem in reference 1. 


3. APPROXIMATE SOLUTION OF THE PROBLEM 


A viscoelastic layer lies on a semi-infinite medium also viscoelastic 
of mass density p:. A force of gravity g per unit mass acts on both 
media and the layer is subject to a horizontal compressive stress (Fig. 1). 
The following assumptions are introduced : 


. 
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(a). The layer behaves like a plate and obeys the simplified equa- 
tions of the plate theory. 

(b). The hydrostatic prestress due to gravity is the only prestress in 
the underlying medium. 

(c). Interfacial adherence of the layer and medium is neglected. 

(d). Both materials are assumed incompressible. 

(e). The weight of the layer has no effect on the folding, i.e., the 
layer density may be assumed equal to zero. 


The justification of assumptions (a) and (b) is established by the 
results in reference 3 which indicate that the factors neglected are not 
too significant for practical purpose. The same remark holds for as- 
sumption (c) on the basis of the analysis carried out in reference 2 by 
taking into account the adherence. Assumption (e) is an approxima- 
tion which for physical reasons may be considered to follow from (a) 
since the plate bending occurs with no variation in thickness. 

It may be expected that the assumption of incompressibility does 
not restrict the results significantly while providing the advantage of 
simplicity in the mathematical treatment. 


4 
YIN Yy 


Yi 


Yy YY Yj iif 


Fic. 1. Compressed layer lying on top of a’ Fic. 2. Coordinates, surface forces, and 
heavy medium. deflections for the half space. 


Consider first the underlying medium as an infinite half-space, the 
surface lying at y = 0, and the y axis being directed downward (Fig. 2). 
Since we neglect any initial horizontal compression P in the underlying 
medium, the equations for the displacement components u and v 
are obtained by putting P = 0 and p = pu, in Eq. 6, hence, 


OS Ov 

0 

OSi2 du 0 

Ox oy Ox 


= = 0. (12) 


Yy 
Since we are dealing with an incompressible material, we may write | 
ou 4 
Ox 
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Introducing this relation in Eq. 11, they become 


0 
(Sir + pipv) + Oy = 0 


(13) 
Siz 
x + ay (Seo + pigv) = 0. 
The boundary conditions (8) at the surface y = 0, when the prestress is 
zero (P = 0), become 

dx 

dx 


(14) 


The operators R and Q of the stress-strain relations (9) defining the 
viscoelastic properties of the underlying medium are designated by R, 
and Q,. For an incompressible material we put (8) R, = ©,e = 0, and 


(15) 
Relations (9) then become 
Sij — = 2Qie;;. (16) 
If we write 
Ss + pigv 
Eqs. 13 are transformed into 


Ox 


Os 12 
Ox 


Relations (16) become 
Sij) — 8558’ = Ques; 
and the boundary condition is 


dF, 
dx 


v= — So 
dx pig 22 


<= 


& 
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OS12 
i (18) 
OSo0' 
S22 
oy 
(19) 

(20) 
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Equations 18 and 19 are those for a medium without initial stress. 
The medium thus behaves as though there were no initial stress except 
for a modified boundary condition 20 which contains an additional 
“buoyancy term’’—pigv. 

The problem of surface deflection of a viscoelastic half-space under a 
sinusoidal vertical load distribution was solved in references 1 and 2. 
For the case of an incompressible material, a vertical load distribution 


aF, 
99°98 lx (21) 


produces a normal deflection v = V cos lx (Fig. 2) related to the load by 
qo = 2Q.1V. (22) 


It was also found that in the absence of a horizontal force at the surface 
(F, = 0) the horizontal displacement at the surface also vanishes for an 
incompressible material. In order to take gravity into account we have 


dF, dF 
just shown that we must replace the boundary force — by om — pigv; 


d. 
hence we must replace qo by go — pigV. Relation 22 then becomes 


go = + pig) V. (23) 


Let us now consider the layer of thickness h. We assume that it bends 
as a viscoelastic plate under an axial compression P and a transverse 
vertical load, 

g = — qcoslx (24) 


equal and opposite in sign to the load applied to the underlying medium. 
We have shown in reference 1 that the plate deflection v satisfies the 
equation 

h* dv dy 


(25) 


where B is the operator 


defining the viscoelastic properties of the layer. For an incompressible 
material R = © and we have 


(27) 


B = 49. | 
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With a sinusoidal deflection v = V cos /x after introducing expressions 
22 and 24 for the load, Eq. 25 becomes 


Phi? + 20.1 = 0. 


(28) 


This is the basic equation of the problem. We note that the mass 
density of the layer does not appear in this equation. This is due 
primarily to our assumption that the deformation of the layer obeys the 
plate bending equation and therefore does not undergo a change of 
thickness. Actually of course a more accurate treatment would show 
an influence of the layer density, but the correction due to this effect 
may be assumed to be small. 

In the next section we shall discuss the significance of Eq. 28 in 
relation to the folding of the layer under the simultaneous action of the 
axial compression P and gravity. 


4. FOLDING OF A VISCOUS LAYER LYING ON A HEAVY VISCOUS MEDIUM 


We shall apply Eq. 28 for an incompressible medium to the par- 
ticular case of a purely viscous layer under a horizontal compression P 
lying on top of another heavy viscous fluid of infinite depth. The 
viscosity coefficient of the layer is denoted by yu, that of the underlying 
medium by wi, and the mass density of the latter by p;. The corre- 
sponding operators are 


Q = up 
29 
Q: 


Relation 28 may then be written 


(30) 


The physical interpretation of this relation lies in the significance of 
the variable p considered as an algebraic quantity. As pointed out in 
our previous work (1, 2, 3) if relation 30 is satisfied, any sinusoidal folding 
of wave length 


L = (31) 


has an amplitude increasing with time ¢, proportionally to the factor 
e?'. We have called the dominant wave length L, that which increases 
at the maximum rate. This dominant wave length corresponds to the 
minimum value of P/up on the right hand side of Eq. 30. This mini- 


4 
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mum depends on two parameters, u/y; and pigh/P. The latter repre- 
sents the influence of gravity, as the ratio of the gravity forces to the 
compressive load P. The relative importance of gravity is therefore 
dependent on the magnitude of the horizontal compression. 

We note that only the density p, of the underlying medium appears 
in the expression for the gravity parameter. As already pointed out 
this is due to our assumption that the layer behaves like a plate in 
flexure and therefore that there is no change in thickness. The domi- 
nant wave length L, of the folding corresponds to values of /h for which 
expression 30 is a minimum. Denoting by /,h this value of Jh, it is 
related to the dominant wave length by the relation 


Lah’ (32) 


20 — 

n 


P 
on 


10 20 


Fic. 3. The combined effect of gravity, Fic. 4. The combined effect of gravity and 
compressive load, and viscosity on the value compressive load on the dominant wave 
of lah. length. 


values of the parameter pigh/P. The case pigh/P = 0 corresponds to 
the absence of gravity. In that case we find the straight line 


lah = (33) 


This is identical with the relation found in reference 1 except for the 
factor 3 in the present formula instead of 6. This is due to the fact that 
the formula of reference 1 applies to a layer embedded in a medium 
on both top and bottom while the present case deals with a layer lying on 


Che value of 1,h is plotted as a function of 4/ —— in Fig. 3 for various 
ignt 
| Lv 
| 
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the surface of a semi-infinite medium. Clearly the latter case amounts 
to the former if we divide by two the viscosity of the embedding medium. 
For pi = 0, that is if the underlying medium is a fluid without viscosity, 
the problem reduces to that of folding of a viscous plate laterally re- 
strained only by the buoyancy of the underlying fluid. This amounts 
to putting uw, = 0 in Eqs. 28 and 30. We derive 


up _ 3 _ 3ogh 1 (34) 


P Fh? P 


The maximum of this expression as a function of /h yields the value 
corresponding to the dominant wave length, that is, 


These are the values plotted along the vertical axis at zero abscissa in 
Fig. 3. From Eq. 35 we may write for the dominant wave length due to 
gravity alone 


La = th pigh’ (36) 


For example if P = gpigh, that is, if the compressive load is nine times 
the weight of a column of the underlying material of height equal to the 
layer thickness we find Ly = 13.4 * h, that is, the dominant wave 
length due to gravity alone is about thirteen times the layer thickness. 

At this point it is of interest to examine the magnitude of the in- 
stability. Following a procedure identical with that used in the pre- 
vious work (2, 3) let us evaluate the factor by which the amplitude of 
the dominant wave length is multiplied for a period of application of 
the compressive load P during which the layer would undergo a com- 
pression of 25 per cent.? This time ¢; satisfies the relation 


Pt, = 4ue = pp (37) 


where the compressive strain is 


e =}. (38) 


We shall examine only the case when the underlying medium has no 
viscosity (4: = 0). The amplification factor is exp(pt,) with 


pigh 


(39) 


3 It is of course understood that the theory is not rigorously applicable to such large de- 
formations and the figures are only given as an indication of orders of magnitude. 
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In Table I we have shown the values of the amplification factor in 
terms of the gravity parameter. 


TABLE I. 
pigh/P exp(pti) 
1 2.11 
1/4 20 
1/7 190 
1/9 850 


It is seen that significant amplification occurs only for P/pigh S$ 7. 
A discussion of the amplification factor in the absence of gravity was 
given in references 2 and 3, and it was found that significant amplifica- 


tion occurred only if u/ui S 70 that is, if Mi a = 0.35. The range of 


the parameters shown in Fig. 3 corresponds to a domain for which the 
amplification factor is significant. 

The combined effect of gravity and viscosity is completely repre- 
sented by the diagram of Fig. 3. We notice that the effect of gravity 
dominates in regions where the wave length due to gravity alone is 
larger than that due to viscosity alone as determined from Eq. 33. For 


example, assuming a viscosity ratio w/p,; = 64, we find Nb = (0.37. 
Mi 

The ratio L,/h of dominant wave length to layer thickness for this case 


is plotted in Fig. 4 as a function of V -t We see that contrary to 
Pi 


the case where gravity is neglected the wave lengths depends on the 
compression P. For large values of P the ratio tends to La/h = 17.2. 
For the smaller values of the compressive load this ratio drops to about 
12, as may be seen from the discussions in references 2 and 3. Below 
this value the magnitude of the instability is such that it loses its 
physical significance, as shown by the values of the amplification factor 
in Table I. 

This means that in order to exhibit appreciable folding a very 
thick layer may require a compressive load beyond the physically 
possible range. ‘This conclusion should be of significance in geological 
applications. 

5. FOLDING OF A LAYER LYING BETWEEN TWO HEAVY MEDIA 

In the previous section we have considered the layer to lie on top 
of another medium of infinite depth. We shall now consider the case 
where another medium of different density and viscosity lies on top of 
the layer. For simplicity the discussion is restricted to materials 
which are purely viscous and incompressible. We designate by p; and 
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wu, the mass density and viscosity of the bottom medium, and by pz» 
and yp, the corresponding quantities for the medium lying on top. The 
layer of viscosity yw and thickness / under a horizontal compression P lies 
between those two media (Fig. 5). 


Fic. 5. Compressed layer lying between two heavy media. 


We may again write Eq. 25 for the layer flexure. However in the 
expression of the total transverse load on the layer g = qocos/x we 
must take into account the restraint due to top and bottom media. 
We therefore write in place of Eq. 23, 


go = [2(Q: + + (p1 — ps)g |V. (40) 
Since the materials are assumed purely viscous, the operators are 
Q: 
= pop. 


The sum Q, + Q:; reflects the fact that the restraint of the two media on 
both sides of the layer is additive while p; — p» corresponds to the fact 
that the effect of gravity is subtractive. Equation 28 for the stability 
is therefore replaced by 


— Phi? + 2(Q: + + (p1 — = 0 


(41) 


(42) 


and Eq. 30 is replaced by 


(43) 


g 


The dominant wave length is determined as before by evaluating the 
value of /h which minimizes this quantity. The equations are formally 
identical with those of the previous section except that the parameter 


h 
ui/mis replaced by and pigh/P is replaced by (pi — 
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We must distinguish two cases. In one case the bottom material is 
heavier than the top one, that is, p: — p2 is positive. The dominant 
wave length is then determined by the same diagrams as in Figs. 3 and 4 
where yu:/u is replaced by (4: + we)/m and pigh/P by (pi — p2)gh/P. 
In the other case the bottom material is lighter than the top one, that is, 
pi — p2 is negative. This case differs essentially from the previous 
one in that an instability due to gravity alone occurs in the absence of 
any compressive stress in the layer.‘ In order to understand this more 
clearly let us write relation 43 in the form 


PPh? + (p2 — pi)gh 


= + + wa) T 


As pointed out above, this value of p gives a folding amplitude propor- 
tional to exp(pt) for the wave length L = 27/l. 

We see immediately that p = , for] = 0, that is, the rate of fold- 
ing increases indefinitely with the wave length. The reader will note 
that we have assumed the deformation to be very small, therefore inertia 
forces have been neglected. Hence, theoretically the dominant wave 
length is infinite. In an actual situation the surrounding medium is not 
infinite and the dominant wave length will be restricted by the thick- 
ness of the surrounding material. This result is of significance in 
geology as it is related to the formation of salt domes. 

There is however the possibility of occurrence of a secondary domi- 
nant wave length due primarily to the compression P. Obviously if 
the effect of gravity disappears, for example if p: = ps, there is a maxi- 
mum value of p for lh equal to 


lsh = (45) 


This relation is found by replacing in Eq. 33 z by (41 + w2)/u. When 


gravity enters into the picture this maximum is displaced. This can 
be shown by writing Eq. 44 in the form 


1/3 
5 = =- 47 
( ) 


* The stability of a stratified heavy viscous fluid has been analyzed from the viewpoint of 
hydrodynamics by S. Chandrasekhar (9) and R. Hide (10). 
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P( ur + 
~ 49 
Gp (49) 


y= 


The problem is to find the maximum value of y considered as a function 
of 6. When we plot y as a function of 6 we must distinguish three 
regions. For 

Kk < 3(4)?? = 0.198, (50) 


the plot shows a maximum and a minimum, as represented by curve 1 
in Fig. 6. For 


Fic. 6. Three typical curves for the dependence of y on 6. 


the maximum and minimum coincide (curve 2) and become a hori- 
zontal inflexion point at 


= V¥} = 0.63. (52) 


For 
> (25) 
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the curve has no maximum or minimum (curve 3). The point M at 
which the maximum of y occurs corresponds to a wave length of maxi- 
mum rate of growth. The values of x and 6 are related at point M by 
the relation 


— 6%) 


267+ 1 


For values of « satisfying the inequality (50) there are two positive 
values of 6 satisfying this relation. The largest value corresponds to 
point M. This largest value of « versus 6 is plotted in Fig. 7. 


| 


Fic. 7. Relation between « and 6 corresponding to point M of Fig. 6. 


From the relationship (Eq. 54) between x and 6 we deduce the value 
of 14h corresponding to the dominant wave length. This is done by 
writing Eqs. 47 and 48 in the form 


6 — pi)gh 
= (55) 


(56) 


+ p2) — pigh 


In these expressions x is a function of 6 through relation (54). They 
define a family of plots for lah versus V3 (us + bo) /m /u with a param- 


eter V(p2 — pi)gh/P which measures the influence of gravity. These 
curves are plotted in Fig. 8, in a way analogous to Fig. 3. We remember 
that the dominant wave length is related to Jah by Eq. 32. The inter- 
rupted line in Fig. 8 corresponds to the case where the maximum van- 
ishes as in curve (2) of Fig. 6, hence to the disappearance of any 
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dominant wave length. The disappearance of the dominant wave 
length occurs for x > 0.198, hence for 


+ ws) — pidgh 


It will be noted that in the case of an inverse density gradient the effect 
of gravity is to increase the dominant wave length. 


02 03 04 
Fic. 8. Combined effect of gravity, compressive load, and viscosity, on the value of Jah for 
the case of an inverse gravity gradient (p2 > p:). 


Finally in connection with the disappearance of the dominant wave 
length if the inequality (57) is satisfied we should remark that this will 
only be true if the surrounding medium is of infinite extent. In ac- 
tuality of course the thickness will have a finite value, and if we neglect 
inertia forces as we have done the dominant wave length caused by the 
inverse gravity gradient will not be infinite but will be determined by 
the thicknesses of the upper and lower media and their own boundary 
conditions. 
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OBLIQUE IMPACT OF A BOUNDED STREAM ON 
A PLANE LAMINA 


BY 
TURGUT SARPKAYA?! 


ABSTRACT 


In this paper, Helmholtz’s free streamline theory is used to determine the con- 
traction coefficients of two jets formed by a finite, two-dimensional stream impinging 
obliquely upon a plane lamina placed between two infinite parallel planes. The 
agreement between the results obtained from experiments and those calculated from 
theory is excellent. Thus, the theory is expected to be useful as a design tool for 
such boundary configurations as those represented by butterfly valves. 


INTRODUCTION 


Many plane irrotational flow patterns which include the formation 
or deflection of free or partly bounded jets can be analyzed satisfactorily 
by the theory of free streamlines devised by Helmholtz (1)? and Kirch- 
hoff (2). The results obtained may be significant not only for the corre- 
sponding two-dimensional flows of real incompressible fluids, but in 
many instances, for their three-dimensional counterparts as well. The 
analysis consists of the definition of successive conformal transforma- 
tions involving a hodograph, or velocity plane, and the application of 
Schwartz (3)—Christoffel (4) transformation. The fact should be em- 
phasized that from the mathematical viewpoint, the direct calculation 
of free-streamline flows is restricted to two-dimensional irrotational 
flows of incompressible fluids which are free from gravitational effects 
and for which the solid boundaries are straight. Hence, the formation 
of two jets by a finite, two-dimensional stream impinging obliquely 
upon a fixed plane lamina placed between two infinite parallel planes 
(Fig. 1) is well suited to free streamline analysis because of the domi- 
nance of inertia and pressure in the establishment of the flow pattern. 

The practical importance of the problem under consideration is 
recognized at once if the similarity of the flow pattern treated herein to 
that encountered in butterfly valves is noted. Butterfly valves in- 
stalled in pipes or penstocks are normally used for one of two basic 
services: either as closure or “‘on-off”’ valves, or for flow rate control or 
“governing.”’ The three major geometric variables affecting flow per- 
formance are the piping arrangement, shape of blade, and the closing 
angle. Available information on contraction and flow coefficients of 
butterfly valves is rather limited and obtained primarily from experi- 
1 Assistant Professor, Department of Engineering Mechanics, University of Nebraska, 


Lincoln, Nebraska. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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ments conducted on two- and three-dimensional models and prototypes. 

This paper includes a demonstration of the effects of geometric 
variables involved and presents the analytical determination of contrac- 
tion and discharge coefficients for a two-dimensional butterfly valve. 
For the present problem since the shear is not a primary factor, for the 
flow accelerates rapidly, the differences between comparable ideal and 
real flows are small. The importance of gravitational effects for the 
coefficient of contraction is already known to be definitely secondary. 
However, the significance of the various restrictions can be shown only 
by comparing the calculated results with those obtained from observa- 
tions of real flows. 


Physical Plane 
Fic. 1. 


STATEMENT OF PROBLEM 


Being fully aware of the restrictions and limitations of the applica- 
bility of this method, the solution of the problem under consideration 
proceeds as follows: If the coefficients of contraction at points B and H 
are C,, and C,2, respectively, as shown in the physical plane (Fig. 1), 
then the main objective of the present investigation becomes the de- 
termination of the relationships: 


C. 


C.. = F,_2(a, B). (1) 


The existing relationships are: 


_¢.2(, _ sine 
Cara sin “) (2) 
C n ve sina (3) 


sin 


Coy = 
B Yj B 
a 
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sin-! — 


in which m and m represent the thickness of the jets at B and H/, respec- 
tively; 7 is the width of channel; 8 is the angle of closure; a@ is the vari- 
able angle; and L is the length of valve plate. The velocities of both 
jets are equal to V,; and the velocity far upstream from the valve plate 
is Vo. 

If z and w are used for the complex variable and the complex poten- 
tial, then the variable ¢ defined as (dw/dz) is given by 


f= -u+w= — g-* (6) 


in which u and v are the velocity components in the x and y directions, 
q is the magnitude of the velocity vector, and @ is the angle of its 


inclination. 
If the variable Q, 


V; 


a=In(-% = (7) 


is introduced, the flow region in the physical plane (Fig. 1) can be 
mapped onto a corresponding region in the Q-plane, asindicated in 
Fig. 2. Either ¢ or 2 must be related to w by means of known methods 
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through still another complex variable designated as ¢t. This relation- 
ship completes the rather complicated functional analysis because, as 
has already been seen, either ¢ or 2 can be expressed in terms of (dw/dz). 

Any polygonal boundary in one complex plane can be transformed 
into the entire real axis of another through application of Schwartz- 
Christoffel theorem, the interior of the polygon being transformed into 
either the upper or the lower half of the second plane. Thus, the out- 
line in the Q-plane is directly transformable into the real axis of the 
t-plane. Therefore, if we open out the Q-plane and lay it along the real 
axis of the ¢-plane as shown in Fig. 3, the interior of the Q-plane trans- 


- +d + 
. . 


E. F H D 


T-Piane 
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forms into the lower half of the ¢-plane. The transforming function 
can be written as follows: 


dt 
ufz — — + N. (8) 


As many factors (¢ — t,) are introduced into the transformation as there 
are vertices of the polygonal boundary in the Q-plane. The value ¢, 
is the locus of the transformed point on the real axis in the ¢-plane, and 
y, is the exterior angle of the polygon at that point in the ¢-plane. 
Three of the t, values can be assigned arbitrarily but with due regard 
to the geometrical characteristics of the Q-plane, and the remainder 
must be given parametric values to be evaluated in terms of given 
quantities, along with M and JN, from the integrated function. 
The desired mapping function between Q and ¢t becomes: 


(t + a)dt 


The transformation equation between the w-plane and the f-plane 
is obtained directly from the complex potentials for sources and sinks, 
point A being considered as a source of strength (27Vo) and points B 
and // as sinks of strengths (2nV,;) and (2mV,), respectively: 


Vo 


n Ln(t — 1) + mLn(t + 1) — + a) (10) 
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_ nV; sin a 


sin B 


X {Ca Ln(t — 1) + Cer Ln(t + 1) — (Cor + Coz) Ln(t + a)} (10a) 


and w is connected to z by the relation 


> etaw. (11) 


The integral given by Eq. 9 can be decomposed into two component 
integrals: 


tdt adt 


EXAMPLES 


Let us consider the transformation of Legendre (5): 


sin’? + bo 
Co sin? @ + do 


(13) 


do, bo, Co, and do can be so chosen that the above integrals simplify and 
reduce to known integrals. It should be noted that these integrals, 
depending on the values of ¢ require separate treatment for each interval 
under consideration. 


Case I 


In the intervals DE, and FE:, that is, for t > d and t < —f, 
one has 


_ @+f)sin? +f) 
(d +f) sin? — (1 +f) 


=sin +S) 0 < 2/2, 
(¢—1)¢@ +f) 


di = — +f)(1 +f) sin cos 
[(d +f) sin? ¢, — (1 +f) 


é 
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If Eqs. 14 and 16 are replaced in Eq. 12, one obtains 


f. 2M + f) sing, +f) 
d) (d +f) sin?¢: — (1 +f)}{ v1 — ky sin? 
2Ma dd. 
= N. (17 


Integral 17 can be separated into a first and third kind elliptic integral : 


~V(1+f)(1+d) 


Il, (¢:, n, k) is the third kind elliptic integral of Legendre, and may be 
written as 


d¢, 
", k) = Sa + n sin? — k? sin? 


The parameter is taken as real, and may take any value, positive or 
negative. For the intervals under consideration, 


<n, < — 1. 


Legendre’s integral is connected with Jacobi’s by the relation 


sn(e1, k 1) 

II = ——~ JI ;(u, 1K’ 20 

u(r 1 1) ~ ki)dn(e1, ki) i( + ) ( ) 

in which «; is written for F(¢:, ki). sn, cn and dn are elliptic functions 

and K’ is the complete elliptic integral of the first kind for the comple- 
mentary modulus. , is given by the relation 


d +f 2 KK! 
and 
L+f d-1 ,, 


= d+f’ cne, = d+f’ dn’e, = ret 


In order to simplify the notations, k1), cn(e:, and dn(e, R:) 
will be denoted by sme,, cne,, and dne,, respectively. 
The solution of Jacobi’s integral is given by 


u,cne,dne, H(¢€, — 1) 


IT at 1K’) = uiZ(€1, + + >Ln 4 (21) 


| 
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in which H/(e€; + u;) is Jacobi’s Theta-function (6), and 


+f) +d)’ 


Hence Legendre’s integral becomes: 


H(e, — 
+ 


II ni, ki) = u,Z (€1, ki) + Ln (22) 


From Eqs. 18 and 22, we have 


SNE, 


—{(1 +a)u - 1) — 


Cne, -dneé, 


H(e, — | 
1 
x ky) + + N (23) 


which is the solution sought for the intervals under consideration. 
The constants M and N are determined from the comparison of Q- and 
t-planes for the corresponding values and were found to be 1 and ta, 
respectively. 


Case II 
For the interval /7B, that is, |t| < 1, if t is taken as 


— 2f sin? ¢. + (1 +f) 
‘2sin? ¢, — (1 +f) 


and the foregoing manipulations repeated, one obtains 


_ 2M 
V¥(1+d)(1+f) 


in which 2, = — -———. which varies between zero and — k.?. Writing 


— ko) = one has 


= F(g2,k 


+= (€2, Rx) + 3Ln us) (26) 


€2 + ue) 


a= 
i 
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x | wale) + 3Ln (27) 


in which 


(1+f)(1 +d) 2(t +f) 


(d — 1)(f — 1) 
(d+ 1)(f +1) 


6(€ + u) also is Jacobi’s Theta-function and connected to H(e + u) by 
the relation 6(z + 1K’) = iBH(z) where 


K and K’ are complete elliptic integrals of the first kind connected by 
the following relationships: 


K = sn-(1,k), = sn—“(1, k’) = sn—(1, V1 — 


Case IIT 
For the interval BD, that is, 1 


<t <d,if tis taken as 


(d — 1) sin? + (1 +d) 


one obtains 


2M 


Repeating the intermediate manipulations, one has 


| J. F. 1. 
and 
| Snes = vj ra’ = = 
eK’ wis : 
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2M 
V(1+f)(1 +d) 


SNE3 
Cn €3° dne; 


{(a +1)u;+2 


x | + 4Ln +N 


€3 + 
in which 


= (d 1) (f 1) = k,’?, = 2(d 


(d + 1)(f +1) +4) 


and 
d—1 
d +1’ 


Cn; = 1 1’ dn; = 


(d + 1)(t — 1) 
(d — 1)(t+ 1) 


= 


Case IV 


For the interval FH, that is, — 1 > ¢ > — f, iftistaken as 


(f — 1) + +1) 


one obtains 


2M 
+a) (d + a)u, — (d+ (+. (32) 


in which 


= 0, = kysn?(te,, k,) = ksn?(e,, ky’) /cn*(€4, k,’) 


ny =ktan?6, sind = sn(e,, ky’). 


Q=u,+- iusZ (ie,) — tes) | 


dn(e, k 4) + tes) 
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From Eqs. 32 and 33, one has 


2M 
Q 
V(1+d)(1 +f) 


sn(€4, ky’) -cn (es, ky’) 


dn(e,, Ry’) 


(a — f)u, — (d +f) 


| +N _ (34) 
in which 


= F (4, Ra), = k,”? = (1 +d) 


k? = (d+1)(f + 1)’ sn k,’) cn(e,, ky’) = 


and 


dn(e, ky’) = 


€, = sn- 


sin (tan- 


Intermediate steps in the solutions of the above integrals have been 
omitted for the sake of brevity. In order to obtain z from Eq. 11, dw 
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should be written as 


nV; _ sin Ce fat 


The rest of the integration could only be performed by numerical 
methods. The best and well-known method for this purpose is King’s 
(7) A. G. M. scale method. — If fora fixed value of 8, numerical values are 
assigned to a, f and d, satisfying the boundary conditions in the trans- 
formed planes, corresponding simultaneous values may be obtained for 
Cea, Ceo and a. 

The results or the functional relationships obtained from these 
equations are shown in Fig. 4. The significance of this figure is obvious. 
The contraction coefficient C.; varies greatly with a, but very little with 
the angle of complete closure 8, whereas the contraction coefficient C-2 
varies appreciably with 8, but very little with a, for the values of a 
greater than approximately 15°. But for small values of the variable 
angle a, C.2: varies extremely fast and approaches unity. The values of 
the corresponding angles (corresponding to a < 10°) used in evaluating 
the elliptic functions approach their limit so closely that difficulties are 
encountered because of limitations in the tables available. 


THEORY VERSUS EXPERIMENTS 


In order to compare the theoretical results with those obtained from 
experiments it is necessary to compute the discharge coefficient from 
theoretically determined contraction coefficients since it is the most 
accustomed and widely used coefficient and more readily determined 
experimentally. The rate of flow and the total head are expressed by 


Q = VA = CoAVgAH (36a) 
AH = — (365) 


in which A is the area of the normal cross-section of the approach 
channel, and Ce the coefficient of discharge. Combining Eqs. 36a and 
366 with Eq. 5, the coefficient of discharge is obtained easily as 


(37) 


sin 


Cg values computed from Eq. 37 for 8 = 75° using the contraction 
coefficients obtained in Fig. 4, are plotted in Fig. 5, and the solid line is 
obtained. In the same figure the experimental data taken from refer- 
ence (8) are also plotted. The agreement between the theoretical and 
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experimental results is fairly good, the theoretical discharge coefficient 
being somewhat smaller. 

A similar computation has been made for 8 = 90°, and the resultant 
theoretical curve with the experimental points taken from reference (9) 
are shown in Fig. 6. The agreement between the theoretical and ex- 
perimental results is excellent. 
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STABILITY OF FLOW 


So far we have dealt with the effects of geometric variables on each 
coefficient of contraction and on discharge coefficient. It is seen that 
the coefficients of contraction are functions of a and B only, and do 
not depend on the magnitude of the rate of flow. 

The second point we propose to discuss is the stability of flow for 
various positions of lamina and in connection with this the explanation 
for the scatter of experimental data for very small values of @ or for 
near-full open positions of lamina. The stability of flow depends on 
the stability of the stagnation point or rather specifically speaking of 
the stagnation line on the valve plate. As a becomes smaller, the stag- 


| | + 
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nation line approaches the lower edge of the valve plate and theoretically 
it reaches the edge when a is zero. But, practically the stagnation 
point or line reaches the sharp end of the plate when a becomes about 
15°. To illustrate this, a curve has been drawn computing the location 
of the stagnation point for 8 = 90° and for various values of a, as shown 
in Fig. 7. The dimensionless parameter \ becomes practically zero for 


a = 20°. 


Coordinates of 
Stagnation Point 
For B-90° 


10° 26 30 40 50 60 70 80 90 


Fic. 7. 


In the case of axially symmetrical flow, the locus of stagnation points 
is curvilinear. Therefore, the stagnation points in the vicinity of valve 
axis are very close to the edge of valve plate, whereas the stagnation 
points in the vicinity of the diameter perpendicular to the valve axis are 
relatively far from the edge. Therefore, it is very likely that in the 
case of axially symmetrical flow for a critical value of a the stagnation 
points relatively near the plate edge jump abruptly to the edge of the 
plate changing the discharge coefficient. In the two-dimensional case 
the entire stagnation line jumps abruptly to the edge of the plate for 
small values of a. Fora critical value of a, however, the stagnation line 
is unstable and jumps to the edge and then back alternately resulting 
in unstable flow. 

It seems to be evident from the foregoing explanation that the 
scatter of the experimental data for small values of a@ is largely due to 
the rapid and abrupt variations in the location of stagnation point in 
that region. Moreover, the geometry of valve plate also has an im- 
portant effect on the variations of aforementioned coefficients. Theo- 
retically, the plate is a thicknessless lamina, but practically it has a 
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finite thickness. Therefore, for a blade with a relatively wide tip, the 
point of separation will be at one or the other edge, depending upon the 
relative blade angle. The sudden slight shift of discharge coefficient 
for some valves is due to this transfer of the separation point. 


SUMMARY 


In this paper, Helmholtz’s free streamline theory has been used to 
determine the geometric characteristics of two jets formed by a finite 
stream impinging obliquely upon a plane lamina placed between two 
infinite parallel planes. Comparison of theory and experiment is made 
and the agreement has been found satisfactory. Characteristics of 
submerged jets, prediction of incipient cavitation behind the lamina, 
and the torque for free-discharge valves will be examined in a subse- 
quent article. 
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THE FRANKLIN INSTITUTE 


MINUTES OF THE STATED MEETING 
February 18, 1959 


The Stated Meeting of The Franklin Institute was heid at 8:15 p.m. in the Lecture Hall. 
Mr. Wynn Laurence LePage, President, called the meeting to order with approximately 300 
members and guests in attendance. 

The President stated that the minutes of the Stated Meeting of December 17, 1958 were 
printed in the January 1959 issue of the JourRNAL. There being no corrections or additions, 
the minutes were approved as published. He announced that the minutes of the Annual 
Meeting of January 21, 1959 will be printed in the February 1959 issue of the JouRNAL and 
will be presented for approval at the next Stated Meeting. 

Mr. LePage announced that all members of The Franklin Institute who had completed 
25 years of continuous membership last year had been invited as our guests this evening. He 
called upon Mr. Morton Gibbons-Neff, a member of the Board of Managers and Chairman of 
the Membership Committee, to present the names of these Silver Members of 1958. Mr. 
Gibbons-Neff read the following names: 


Charles A. Bareuther Annetta R. Masland 
George E. Crofoot Ellen-Jarden Nolde 
George S. Gardner J. Livingston Poultney 
Walter H. Hartung George A. Smith 

Karl S. Howard Floyd T. Tyson 
Edward R. Hughes Wilfred H. Wickersham 


Mr. LePage then announced that this evening would mark the Fourth Annual Philip C. 
Staples Lecture. This lecture was established by friends of Mr. Staples, late president of 
The Franklin Institute and the Bell Telephone Company, as a most fitting memorial to a man 
who had done so much not only for The Franklin Institute but for the community in which 


he lived. 

He introduced the speaker of the evening, Mr. H. Thomas Hallowell, Jr., President of 
Standard Pressed Steel Company of Jenkintown, Pennsylvania, whose subject, ‘Now It’s Got 
To Work” explained why high reliability design and execution is so necessary. Mr. Hallowell’s 
talk placed emphasis on the reliability factor, so necessary today because of the manifold parts 
and units that go into making up the equipment used in the home, in industry, in transporta- 
tion, in our military operations and in other facets of our every-day life. An interesting film 
produced by Standard Pressed Steel Company illustrated his talk. Mr. Hallowell spoke of 
the desirability of The Franklin Institute, as a neutral organization, taking the lead in ar- 
ranging for a Symposium on this subject, to include Government agencies, industries and 
technical societies. 

Interest in Mr. Hallowell’s subject was very evident by the half-hour question period and 
the applause following his talk. After thanking Mr. Hallowell for his very interesting talk, 
the President adjourned the meeting at 9:55 p.m. 


WILL1aM F. Jackson, JR. 
Secretary 
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MEMBERSHIP 
1958 SILVER MEMBERS HONORED 


At the Stated Meeting on February 18, 1959, the Institute honored twelve of its members 
who had completed 25 consecutive years of membership. Prior to the meeting, the four new 
Silver Members who were able to attend were guests of honor at a dinner and reception. 
During dinner, Wynn Laurence LePage, President of the Institute, introduced these honored 
guests: George S. Gardner, Karl S. Howard, Miss Annetta R. Masland and Floyd T. Tyson. 
As souvenirs of the occasion, the new Silver Members received initialled paper weights con- 
taining 1958 quarters, presented to them by Morton Gibbons-Neff, Chairman of the Mem- 
bership Committee. Mr. Neff read the names of the twelve new Silver Members at the Stated 
Meeting following the dinner. 


BRIEF BIOGRAPHICAL NOTES ON THE 1958 SILVER MEMBERS 


Cuar_es A. BAREUTHER—Chairman of the Physics Department of the Central High 
School of Philadelphia and Professor Emeritus of the Evening College of Drexel Institute 


of Technology. 

GerorGE S. Croroot—Retired Professor of Mechanical Engineering at the University of 
Pennsylvania. 

GEORGE S. GARDNER—Research Chemist at Amchem Products Inc., who has contributed 
articles on corrosion and inhibition to this JOURNAL. 

Water H. HartunG—Professor and Chairman of Chemistry and Pharmaceutical Chem- 
istry at the Medical College of Virginia. 


Kari S. Howarp—Retired Vice President of General Steel Castings Corporation; now 
Consultant to the U. S. Army Ordnance Tank-Automotive Command in Detroit. 


Epwarp R. HuGHes—Electrical Engineer who specializes in patent work ranging from 
refrigerators to weeping dolls. 


ANNETTA R. MASLAND—Retired biology teacher who taught in West Philadelphia High 
School and Overbrook High School. 


ELLEN-JARDEN No_tpeE—Psychologist connected with testing programs in various Phila- 
delphia schools; United Nations correspondent for weekly newspapers. 


J. LivinGstone Pou_ttNey—Retired Consulting Engineer for the Brill Corporation. 
GerorGeE A. SmirH—Patent lawyer and Partner in the firm of Busser, Smith and Harding. 
FLoyp T. Tyson—Professor of Chemistry at Temple University. 


Witrrep H. WicKERSHAM—Resident of Nutley, N. J. 
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The following elections were approved on February 5, 1959 by the Membership Committee 
by authority given to it by the Board of Managers. 


Jacques D. Bambling 
Francis J. Bonner 
Carl Richard Borson 
Newman Bortnick 
Ralph D. Bradway 
Daniel B. Campbell 
James B. Chambers 
F. W. Chapman, Jr. 
Charles W. Charny 
Ellis E. Comfort 
John A. Cost 

Jack Dickstein 
Hans O. Dietze 


Joseph N. DuBarry, IV 


Kenneth V. Duncan 
William Ehrich 
John L. Frost 

Loren E. Gaither 

H. Paul Gant 
Augustine M. Garcia 


John F. Harvey 


Julian M. Avery 


Seth Bunker Capp ‘11 
Elizabeth B. Chew '36 
J. S. Curley '36 
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NEW MEMBERS 


ANNUAL 


Miro A. Grottger 
John E. Jednacz 
Charles H. Justice 
Nelson J. Katz 
Lucien Katzenberg, Jr. 
Jerome W. Kaufman 
James H. Keegan 
Clyde M. Keitt 
Louis W. Kellman 
Calvin F. Kendig 
Walter S. Kenworthy 


Charles R. King 


Lester Kiriloff 

L. R. Kirkwood 

T. S. Klingeman 
Howard S. Kneedler, Jr 
Karl F. Koehler 

Louis Kolodin 

Charles E. Kunz 

James A. Mallon 

John Massier 


EDUCATOR 
Fred M. Hofkin 
Howard B. Tansley 
NON-RESIDENT 


B. Ward Deutschman 


NECROLOGY 


Edwin H. Dressel '45 
E. Monroe Fisher '52 
Margaret Horstmann Packard '35 


J. Grant McCabe, III 
Richard S. Lee 
Richard S. Loveland 
Thomas J. McCann 
Carl V. McCarty 
Harry Luecke 
William Metcalf 
Howard R. Milligan 
Charles P. Neidig 
Arnold N. Parris 
Charles L. Register 
George F. Rickert 
Henry L. Rothman 
Paxton G. Seabury 
John L. Tetley 
Richard T. Urffer 
Maria Giulia Vitelli 
Irwin Frank Wertleib 
Joseph T. Wolpert, Sr. 
Ernest R. Wortmann 


Hsiu Huo Yuen 


Benjamin A. Oxnard 


Mrs. Benjamin Rush, Sr. '36 


Charles S. Redding '26 
Charles H. Schaefer '43 
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The following are Silver Members, having belonged continuously for twenty-five years 
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or more. The dates are years of election to membership. 


1910 
Richard P. Brown 
Joseph S. Hepburn 
Joseph W. Lippincott 
J. D. Shattuck 

J. M. Weiss 


1912 

James Barnes 

Morris Llewellyn Cooke 
Edwin Elliot 

Walter Palmer 

Julian S. Simsohn 
Howard S. Worrell 


1913 

Sterling H. Bunnell 
Irenee duPont 

S. B. Eckert 

Ralph E. Flanders 

N. E. Funk 

Lionel F. Levy 

E. Mallinckrodt, Jr. 
Frank M. Masters 
Franklin P. McConnell 
William Maul Measey 
W. Chattin Wetherill 
D. Robert Yarnall 


1914 

Russell L. Brinton 
George S. Crampton 
Ernest L. Huff 

K. G. Mackenzie 
George Wharton Pepper 


1915 

Thomas D. Cope 
John C. Cornelius, Jr. 
H. Jermain Creighton 
Thomas W. Elkinton 
Richard Howson 

W. Wallace McKaig 


1916 

Henry B. Allen 
James G. Detwiler 
Zay Jeffries 

Frank H. Sauer 

B. E. Shackelford 
Arthur Synnestvedt 
Charles N. Weyl 


1917 

A. D. Chambers 
Francis J. Chesterman 
Walter S. Crowell 


William J. Fitzmaurice, Jr. 


Rolfe E. Glover, Jr. 
Harold Goodwin, Jr. 
Arthur W. Lowe 
Frank S. MacGregor 
Francis F. Milne, Jr. 
C. H. Quinn 

Sylvan D. Rolle 
Robert L. Wood 
Henry Woodhouse 


1918 
Charles H. Colvin 
Walter O. Snelling 


1919 

Charles E. Brinley 
Henry Colvin, 2nd 

E. A. Eckhardt 
George W. Furness 
Frank H. Griffin 
Hiram S. Lukens 

J. Howard Pew 
Edmund G. Robinson 
L. K. Sillcox 


1920 

Robert Cameron Colwell 
Chester Lichtenberg 

P. S. Lyon 

Haviland H. Platt 

Carl D. Pratt 


1921 
Roland L. Andreau 
Leonard T. Beale 

J. Ed. Brewer 
Arthur L. Day 
Frederic Palmer, Jr. 
Howard Stoertz 
James Stokley 
Edward R. Weidlein 
S. Weinberg 


1922 

C. Douglas Galloway 
Joseph F. Greene 

C. E. K. Mees 


1923 

Henry Clay Gibson 
Howard S. Levy 
William C. Melcher, Jr. 
John F. Metten 
Richard G. Sagebeer 
H. Birchard Taylor 
Francis B. Vogdes 
Joseph J. Vogdes 


1924 

Peter Abrams 
William G. Ellis 

T. R. Harrison 
Valentine Hiergesell 
Karl F. Oerlein 
Edward B. Patterson 
Nicholas G. Roosevelt 
Alfred O. Tate 

James Lloyd Weatherwax 
Alexander Wilson, 3rd 
William Zimmermann 


1925 

Jack G. Binswanger 
Marion Eppley 
John Graham Foley 
J. V. Giesler 

C. R. Kraus 

H. F. Porter 

W. F. G. Swann 


1926 

Charles B. Bazzoni 
Neils Bohr 

O. E. Buckley 
William Findlay Downs 
Henry C. Evans 
Edward L. Forstall 
Lewis H. Hendrixson 
Charles S. Leopold 
Alfred L. Loomis 

J. H. Manning 
Clement B. Newbold 


O. M. Patton 
Dexter N. Shaw 
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1926 (cont.) 
Sydney L. Wright 
Charles S. Wurtz, Jr. 


1927 

T. K. Cleveland 

William DeKrafft 

Hiram B. Ely 

Leslie Griscom 

Alfred Iddles 

John Van Gasken Postles 
Philip Sporn 

Charles A. Stanwick 
Charles A. Young 


1928 

William R. Dohan 
Walton Forstall, Jr. 
C. W. Hansell 
Walter A. MacNair 
Fred P. Nabenhauer 
Joseph N. Pew, Jr. 
William Rogers, Jr. 
Frank N. Speller 


1929 

N. E. Bonn 

Bernard Davis 

Rupen Eksergian 
Edward C. Haines 
Malcolm Lloyd Hayward 
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1929 (cont.) 

Randal Morgan, III 
Harold F. Pitcairn 
Lawrence Saunders 

A. Kelvin Smith 

I. Melville Stein 
William Stericker 
Harden Franklin Taylor 
Jacques L. Vauclain 
Albert J. Williams, Jr. 
Arthur M. Wilson, Jr. 
C. T. R. Wilson 


Enos E. Witmer 


1930 

Mildred Allen 
William M. Barret 
Aldus Fogelsanger 
John Gardiner, Jr. 
Morton Gibbons-Neff 
Norman R. Gibson 
Carroll Hodge 
Walter J. Lehman 

R. P. Long 

S. B. Morehouse 
William C. Peterman 


G. E. Wiltbank 


1931 
Charles E. Berger 
Carl Breer 


1931 (cont.) 

Daniel B. Curll, Jr. 
Edward Dawson 

J. Bennett Hill 
Frederick C. Holtz 
Penrose R. Hoopes 
William Melas 
Parry H. Moon 
Julius Neumueller 
Laurence P. Sharples 


1932 
Arthur Fleischer 
William A. Heine 


Wynn Laurence LePage 


Robert McLean 
E. Burke Wilford 


1933 

Charles A. Bareuther 
George E. Crofoot 
George S. Gardner 
Walter H. Hartung 
Karl S. Howard 
Edward R. Hughes 
Annetta R. Masland 
Ellen-Jarden Nolde 


J. Livingstone Poultney 


George A. Smith 
Floyd T. Tyson 


Linn Helander Abraham D. Caesar Wilfred H. Wickersham 


The following are Gold Card Members, having belonged continuously for fifty years or 
more: 


Frank Shaw Clark '07 
W. M. C. Kimber '08 
Theobald F. Clark '09 
Walton Clark, Sr. '09 
James McGowan, Jr. '09 


Charles H. Howson '03 
Frank S. Busser '04 

C. Mahlon Kline '04 
George R. Hall '05 
Darthela Clark '07 


Fred H. Colvin ’88 
G. H. Clamer '91 

J. H. Granberry '95 
Edward Woolman '95 
E. A. Muller '99 


James M. Caird '01 


T. Edward Ross '09 
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LIBRARY 


The Committee on Library desires to add to the collections any technical works that 
members may wish to contribute. Contributions will be gratefully acknowledged and placed 
in the Library. Duplicates received will be transferred to other libraries as gifts of the donor. 
365 items have been added the past month. 

Photostat service. Photostat prints of any material in the collections can be supplied 
on request. 

The Library and reading room are open from 9 A.M. until 5 p.m. on Mondays, Tuesdays, 
Thursdays, Fridays and Saturdays; and 2 P.M. until 10 p.m. on Wednesdays. 


RECENT ADDITIONS 
AERONAUTICS 


Pore, ALAN. Aerodynamics of Supersonic Flight. Ed. 2. 1958. 
Zucrow, Maurice J. Aircraft and Missile Propulsion. Vols. 1-2. 1958. 


ARCHITECTURE AND BUILDING 


PARKIN, PETER HUBERT AND HUMPHREYS, HENRY ROBERT. Acoustics, Noise and Buildings. 
1958. 

SEVERNS, WILLIAM HARRISON AND FELLOWS, JULIAN ROBERT. Air Conditioning and Re- 
frigeration. 1958. 

ASTRONOMY 

BaiLy, Francis. Astronomical Tables and Formulae Together with a Variety of Problems 
Explanatory of their Use and Application. 1827. 

FERGUSON, JAMES. Astronomy Explained upon Sir Isaac Newton's Principles. Ed.9. 1794. 

HERSCHEL, JOHN FREDERICK WILLIAM. Outlines of Astronomy. 1876. 

HuMBoLptT, FriepricH HEINRICH ALEXANDER VON. Cosmos; a Sketch of a Physical De- 
scription of the Universe. Vols. 1-4. 1850-1852. 

KLINKERFUES, ERNST FRIEDRICH WILHELM. Theoretische Astronomie. 1871. 

LAPLACE, PIERRE Simon, MARQuIS DE. The System of the World. Vols. 1-2. 1830. 

Ler, Tuomas J. A Collection of Tables and Formulae Useful in Geodesy and Practical 
Astronomy, Including Elements for the Projection of Maps. 1849. 

NARRIEN, JOHN. Practical Astronomy and Geodesy. 1845. 

NEwcos, SIMON AND HOLDEN, EDWARD SINGLETON. Astronomy for Schools and Colleges. 
1879. 

Denison. A Compendium of Astronomy. Ed. 2. 1842. 

PONTECOULANT, PHILIPPE GUSTAVE DouLcEeT, GRAF DE. Théorie Analytique du Systemé 
du Monde. Ed. 2. Vol. 1. 1856. 

PoNnTIFICAL ACADEMY OF SCIENCE. Stellar Populations. 1958. 

SCHWARZCHILD, MARTIN. Structure and Evolution of the Stars. 1958. 

STOCKWELL, JoHN N. Theory of the Moon's Motion. 1881. 

Wuiston, WiLL1AM. Astronomical Lectures. Ed. 2. 1728. 


AUTOMATIC CONTROL 
CosGrirF, RoBert LizeN. Nonlinear Control Systems. 1958. 
EcKMAN, Donatp P. Automatic Process Control. 1958. 
Hotzpock, WERNER G. Automatic Control, Principles and Practice. 1958. 
Jury, Ipranim. Sampled-Data Control Systems. 1958. 
SUTHERLAND, Ropert L. Engineering Systems Analysis. 1958. 


BIOGRAPHY 


BLOCHMAN, LAWRENCE GOLDTREE. Doctor Squibb; the Life and Times of a Rugged Idealist. 
1958. 
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BIOLOGICAL CHEMISTRY 


BriopHysicaL Society. ist SyMpostuM, CAMBRIDGE, MaAss., 1958. Microsomal Particles 
and Protein Synthesis. 1958. 


BOTANY 
CocHRANE, VINCENT W. Physiology of Fungi. 1958. 


CHEMISTRY AND CHEMICAL TECHNOLOGY 


BAGNALL, K. W. Chemistry of the Rare Radioelements. 1957. 

Battista, O. A. Fundamentals of High Polymers. 1958. 

BLUMENTHAL, WARREN B. The Chemical Behavior of Zirconium. 1958. 

Boxitz, Davip F., Ep. Colorimetric Determination of Nonmetals. 1958. 

Bonnar, RoBert U.; DimBat, MARTIN AND StROsS, FRED H. Number-Average Molecular 
Weights. 1958. 

Bovey, FRANK A. The Effects of Ionizing Radiation on Natural and Synthetic High Poly- 
mers. 1958. 

Cook, RoBert P., Ep. Cholesterol. 1958. 

Gas Chromatography. 1958. 

GorRuP VON BESANEZ, EUGEN FRANZ. Lehrbuch der Chemie fiir den Unterricht auf Univer- 
sitaten und mit Besonderer Beriicksichtigung des Standpunktes Studirender Mediciner. 
Vol. 1-2. 1859-60. 

GunstTonE, F. D. An Introduction to the Chemistry of Fats and Fatty Acids. 1958. 

MOELLER, THERALD. Qualitative Analysis. Ed. 1. 1958. 

NEILANDS, J. B. AND Stumpr, Paut K. Outlines of Enzyme Chemistry. Ed. 2. 1958. 

PENNER, S.S. Chemistry Problems in Jet Propulsion. 1957. 

E. EMMet. Organic Chemistry of Bivalent Sulfur. Vol. 1. 1958. 

SCHWARZENBACH, GEROLD. Die Komplexometrische Titration. Ed. 3. 1957. 

Society oF CHEMICAL INDUSTRY. Industrial Carbon and Graphite. 1958. 


COMPUTERS 


ALT, Franz L. Electronic Digital Computers. 1958. 
Karp_us, WALTER J. Analog Simulation. 1958. 
NEUMANN, JOHANN LupwIiG von. The Computer and the Brain. 1958. 


ELECTRICITY AND ELECTRICAL ENGINEERING 


CAvUER, WILHELM. Synthesis of Linear Communication Networks. Ed. 2. 1958. 

Dummer, G. W. A. Fixed Resistors. 1956. 

HARRINGTON, ROGER F. Introduction to Electromagnetic Engineering. 1958. 

HIGONNET, RENE A. AND GREA, RENE A. Logical Design of Electrical Circuits. 1958. 

RapIo-ELECTRONICS-TELEVISION MANUFACTURERS ASSOCIATION. Proceedings of the 2nd 
RETMA Conference on Reliable Electrical Connections, 1956. 1957. 


ENGINEERING 


Boyer, WALTER C. AND ABRAMS, JOEL I. Influence Lines for Continuous Beams. 1958. 

DuRELLI, AuGusto J.; Puitirps, E. A. Anp Tsao, C. H. Introduction to the Theoretical 
and Experimental Analysis of Stress and Strain. 1958. 

HAMMOND, Ro_t. Water Power Engineering and Some Electrical Problems. 1958. 

Recent Advances in the Engineering Sciences. 1958. 


FOOD 


ANDERSON, Oscar Epwarp. The Health of a Nation; Harvey W. Wiley and the Fight for 
Pure Food. 1958. 
Jacoss, Morris Boris. The Chemical Analysis of Foods and Food Products. Ed. 3. 1958. 
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GEOLOGY 


Busk, Henry Goutp. Earth Flexures. 1957. 
GuiLcHerR, AnpR&. Coastal and Submarine Morphology. 1958. 
HernricH, EBERHARDT WILLIAM. Mineralogy and Geology of Radioactive Raw Materials. 


1958. 
Mason, Brian. Principles of Geochemistry. Ed. 2. 1958. 


HOROLOGY 


AprRIEN, E. Aux Horlogers. 1888. 

BerLLarp, ALFRED. La Montre. 1907. 

Erizzo, NicoLco. Relazione Storico-Critica della Torre dell 'Orologio di S. Marco in Venezia. 
1860. 

Foucuer, M. P. Manuel d’Horlogerie Eleve des Fabriques de Geneve. 1850. 

HEssEN, ANDRE. Memoire sur l'Horlogerie. 1784. 

HILLMANN, Bruno. Die Reparatur Komplizierter Taschenuhren. 1911. 

Lioyp, H. ALAN. Some Outstanding Clocks over Seven Hundred Years 1250-1950. 1958. 

SirE, GeorGes. L’Horlogerie a l'Exposition Universelle de 1867 a Paris. 1870. 


INDUSTRIAL MANAGEMENT 
BARNES, RALPH Mosser. Motion and Time Study. Ed. 4. 1958. 
MANUFACTURE 


Apps, E. A. Printing Ink Technology. 1958. 
CoLeman, D.C. The British Paper Industry, 1495-1860. 1958. 
HERRMANN, ERHARD. Handbuch des Stranggiessens. 1958. 


MATHEMATICS 


ApLER, IRviING. The New Mathematics. 1958. 

Court, NATHAN ALTSHILLER. Mathematics in Fun and in Earnest. 
Cox, D. R. Planning Experiments. 1958. 

FrerGusON, ROBERT O. AND SARGENT, LAUREN F. Linear Programming. 1958. 

ForsyTHE, GEORGE ELMER AND ROSENBLOOM, PauL C. Numerical Analysis and Partial 


1958. 


Differential Equations. 1958. 

FRANKLIN, Pattie. An Introduction to Fourier Methods and the Laplace Transformation. 
1958. 

Hurewicz, WitoLp. Lectures on Ordinary Differential Equations. 1958. 

KAPLAN, WILFRED. Ordinary Differential Equations. 1958. 

Onoe, Morro. Tables of Modified Quotients of Bessel Functions of the First Kind for Real 
and Imaginary Arguments. 1958. 

PaGe, CHEsTER HALL. Algebra of Electronics. 1958. 

Ricumonp, A. E. Calculus for Electronics. 1958. 

SOKOLNIKOFF, IVAN STEPHEN AND REDHEFFER, R. M. Mathematics of Physics and Modern 
Engineering. 1958. 

Some Aspects of Analysis and Probability. 1958. 

TayLor, AnGus Introduction to Functional Analysis. 1958. 

WaLLace, ANDREW H. Homology Theory on Algebraic Varieties. 1958. 


MECHANICAL ENGINEERING 


Heat Exchangers; Applications to Gas Turbines. 1958. 


HRYNISZAK, W. 
Proceedings of the Conference on Lubrication 


INSTITUTION OF MECHANICAL ENGINEERS. 


and Wear. 1958. 
METALLURGY 


CONFERENCE ON THE PHysICAL CHEMISTRY OF IRON AND STEELMAKING. The Physical 
Chemistry of Steelmaking; Proceedings. 1958, 
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The Metal Molybdenum. 1958. 
New Jersey Zinc Co. Zinc Oxide Rediscovered. 1957. 


METEOROLOGY 


Bates, Davip Ropert. Earth and its Atmosphere. 1958. 
DEFANT, ALBERT. Ebb and Flow; the Tides of Earth, Air and Water. 1958. 


MINING ENGINEERING 


Progress in Mineral Dressing. 1958. 


MISSILES AND ROCKETS 


BOUpELER, WERNER. Operation Vanguard. 1958. 

Dow, RicHarp B. Fundamentals of Advanced Missiles. 1958. 

HaAvLey, ANDREW GALLAGHER. Rocketry and Space Exploration. 1958. 

SPECIAL COMMITTEE OF THE INTERNATIONAL GEOPHYSICAL YEAR. Annals of the Inter- 
national Geophysical Year. 1958. 

U.S. Arr Force. The United States Air Force Report on the Ballistic Missile. 1958. 


NUCLEAR ENGINEERING 
Elsevier's Dictionary of Nuclear Science and Technology, in Six Languages: English/American, 
French, Spanish, Italian, Dutch, and German. 1958. 
Hausner, Henry H. Powder Metallurgy in Nuclear Engineering. 1958. 
STEPHENSON, RICHARD. Introduction to Nuclear Engineering. Ed. 2. 1958. 


OPTICS 


Bott, R. H. AND Otners. ‘Tables of Light-Scattering Functions Relative Indices of Less 
than Unity and Infinity. 1958. 

Gray, PETER. Handbook of Basic Microtechnique. Ed. 2. 1958. 

Have, Artuur J. Interference Microscope in Biological Research. 1958. 

SMITH, SIGMUND L. Ore Microscopy. Ed. 2. 1957. 


PETROLEUM TECHNOLOGY 


Brooks, BENJAMIN TALBOTT AND OTHERS, ED. The Chemistry of Petroleum Hydrocarbons. 
1954-55. 
PHOTOGRAPHY 
ALLEN, Roy Morris. Photomicrography. Ed 2. 1958. 
Baines, H. The Science of Photography. 1958. 
SrpLey, Louts WALTON. Photomechanical Halftone. 1958. 


PHYSICS 


CutsHoum, J. S. R. AND Borne, A. H. pe. An Introduction to Statistical Mechanics. 1958. 

DunGey, J. W. Cosmic Electrodynamics. 1958. 

ELFron, ALEXANDER. Mechanics. 1958. 

FROHLICH, HERBERT. ‘Theory of Dielectrics; Dielectric Constant and Dielectric Loss. Ed. 2. 
1958. 

Jackson, J. D. The Physics of Elementary Particles. 1958. 

LAPLACE, PrERRE Simon, MaARQuis DE. A Treatise upon Analytical Mechanics. 1814. 

Nuclear Masses and their Determination. 1957. 

SINNOTT, MAuRICE J. The Solid State for Engineers. 1958. 

STEPHENSON, G. AND KILMIsTER, C. W. Special Relativity for Physicists. 1958. 

TempLe, G. An Introduction to Fluid Dynamics. 1958. 

THIRRING, WALTER E. Principles of Quantum Electrodynamics. 1958. 
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PROGRAMMING 
Gass, Saut |. Linear Programming: Methods and Applications. 1958. 
METZGER, Ropert W. Elementary Mathematical Programming. 1958 
Vazsony1, ANDREW. Scientific Programming in Business and Industry. 1958. 


SANITARY ENGINEERING 


Baspitt, HAROLD E. AND BAUMANN, E. RoBERT. Sewerage and Sewage Treatment. Ed. 8. 
1958. 

McCaBE, JOSEPH AND ECKENFELDER, W. W., ED. Biological Treatment of Sewage and In- 

Vol. 1-2. 1956-58. 


dustrial Wastes. 


SCIENCE 


HANsoN, Norwoop RussELL. Patterns of Discovery. 1958. 
ScHWARTZ, GEORGE I. AND BisHop, Pattie W., ED. Moments of Discovery. 1958. 


SCIENTIFIC INSTRUMENTS 
ADLER, IRvING. ‘The Tools of Science; from Yardstick to Cyclotron. 1958. 


STATISTICS 


GuMBEL, Emit Jurus. Statistics of Extremes. 1958. 
HAGGARD, Ernest A. Intraclass Correlation with the Analysis of Variance. 1958. 


wooD 
Holztechnologisches Handbuch. 1949-58. 


VORREITER, LEOPOLD. 


DO YOU KNOW OF THE FRANKLIN INSTITUTE’S 
NEED FOR FINANCIAL SUPPORT 


Since 1824 The Franklin Institute, through wise and conservative management, 
has been partially self-supporting. 

As an old institution, established in the public interest, The Franklin Institute 
has tried to earn most of its cost by diligent and productive work. Like many 
other institutions of a non-profit nature, The Institute has for many years found it 
increasingly difficult to perform the broad educational services of which the growing 
nation and territories are ever in need 

The educational service rendered by The Institute begins with the early grades 
of our elementary schools and continues throughout an individual's professional and 
industrial career. This period represents a span of fifty years. 

Despite an erroneous public impression gained, perhaps, by the excellence and 
wide range of services performed, that it is a wealthy, richly endowed organization, 
The Institute always needs additional funds to help it perform its broad educational 
program. 

There is a warm satisfaction in giving financial support to an organization that 
will use it constructively to help advance and improve our type of civilization. 

Your gift or bequest, large or small, will be deeply appreciated. We shall be 
pleased to give any additional information regarding gifts or memorials. 

Where property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

Will you please write to: The Secretary, The Franklin Institute, Benjamin 
Franklin Parkway at Twentieth Street, Philadelphia 3, Pennsylvania, for any addi- 
tional information. 
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BOOK REVIEWS 


MATHEMATICAL THEORY OF COMPRESSIBLE 
FLurp FLow, by Richard von Mises, com- 
pleted by Hilda Geiringer and G. S. S. 
Ludford. 514 pages, diagrams, 6 X 9 in. 
New York, Academic Press Inc., 1958. 
Price, $15.00. 


It is no mean task to combine rigorous anal- 
ysis with easy readability in presenting the 
mathematical theory of compressible fluid flow. 
This has been successfully achieved in this 
book. The first three chapters, written by the 
late Professor Richard von Mises before his 
death, provide the keynote for the rest of the 
book which was completed by his wife, Hilda 
Geiringer, together with G. S. S. Ludford. 

In Chapter I, the basic differential equa- 
tions of fluid flow in three dimensions are care- 
fully derived in a completely general manner, 
including the effects of viscosity and heat con- 
duction. The origin and necessity of each 


equation is pointed out, together with its im- 
plications, thus laying a solid foundation for 
the superstructure of the mathematical theory 


of compressible fluid flow. The underlying 
assumptions are examined, as for example, 
that of a perfect gas, and in addition it is 
shown how a non-perfect gas can be treated 
analytically. The author uses the novel con- 
cept of a specifying equation which clarifies 
the role of thermodynamics in the mechanics 
of compressible fluids. 

The mathematical analysis is thorough, yet 
simply presented, and no intervening steps in 
any derivation are omitted on the assumption 
that they are obvious to the reader. The 
author carefully paves the way for the intro- 
duction of vector notation and the use of in- 
tegral theorems in a manner that will be ap- 
preciated by the student who is barely ac- 
quainted with these subjects. A typical ex- 
ample will be found in Article 6.1. Here the 
author derives Stokes theorem in vector nota- 
tion before introducing the Vortex Theory of 
Helmholtz and Kelvin. Of all the different 
derivations of Stokes theorem that this re- 
viewer has seen, not one surpasses this for 
clarity and simplicity. 


Small perturbation theory, which plays an 
important role in aerodynamics is not con- 
sidered further in this book after Article 4 in 
the first chapter. It was the author’s inten- 
tion to discuss extensively the approximation 
theories in the second part of his work; this, 
unfortunately, he did not live to do. 

The general theorems of vorticity and irro- 
tational motion are presented in Chapter II, 
along with the relationships of steady flow. 
Two noteworthy features of this chapter are 
the thorough, yet lucid, treatment of the hodo- 
graph representations and the method of 
charicteristics. 

The third chapter specializes the general 
equations previously described to those of one 
dimensional flow. Solutions are given for flow 
with viscosity and heat conduction and of the 
non-steady flow of an ideal fluid. Examples 
of simple waves and the theory of shock phe- 
nomena bring this chapter to a close, and with 
it, the work of Richard von Mises. 

The last two chapters were written by the 
co-authors of the book. The further develop- 
ment of the theory of plane steady potential 
flow is to be found in Chapter IV, which was 
written by Hilda Geiringer, who also wrote 
the text on the development of Chaplygin’s 
Method and the article on Transonic Flow in 
Chapter V. The bulk of Chapter V is an ex- 
cellent treatment of shock theory and was 
written by G. S. S. Ludford. The clarity of 
presentation combined with thoroughness of 
analysis which characterizes the first three 
chapters is also present in the last two. Asa 
matter of fact, the smoothness of style uni- 
formly present throughout the book leads one 
to suspect that the co-authors played a large 
part in achieving this. 

The book was written as an advanced text- 
book for both graduate students and research 
workers. Its complete lack of exercises for 
the student prevent one from prescribing it as 
a classroom text, but its other qualities make 
it indispensable to the student who wishes to 
gain deeper insight into the mathematical 
theory of compressible flow. The research 
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worker will find this book a valuable addition 
tohis library. The practising gas dynamicists 
who had their formal training in engineering 
schools will not find themselves overwhelmed 
by the mathematics in this book; those who 
entered the field of gas dynamics after com- 
pleting their education in physics or mathe- 
matics will find all the mathematical exacti- 
tude they desire. 
ARTHUR A, Ezra 
The Martin Company 


VacuuM METALLURGY, edited by R. F. Bun- 
shah. 472 pages, illustrations, 6 x 9 in. 
New York, Reinhold Publishing Corp., 
1958. Price, $12.50. 


This book consists of a series of 28 lectures 
delivered at the Department of Metallurgical 
Engineering of the New York University, Col- 
lege of Engineering in June, 1957. It is to be 
the first in a new series on Materials Tech- 
nology. The book comes at a time when 
vacuum metallurgy is of considerabie interest 
in the solution of problems involving pure 
metals and alloys. There have been great 
changes in this field in the last ten years. Un- 
fortunately, there are still a great number of 
researchers who have not followed the tech- 
nology and who still use vacuum techniques 
that they used in school. This book should 
be of great value to bring them, as well as 
others, up to date in this fast moving field. 

It is difficult to collect the individual lec- 
tures of a group of men into one text and still 
maintain a balanced coordinated coverage of 
Dr. Bunshah has done a credit- 
able job in this respect. The format is logical 
and the coverage excellent. The 28 lectures 
are proportioned into ten parts dealing with 
vacuum equipment, thermodynamics and ki- 
netics, arc melting, induction melting, electron 
bombardment melting, degassing, distillation, 
applications, analytical techniques, and future 
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trends. 

The section on applications occupies about 
a fifth of the text. It covers current practice 
in specified areas and also gives comparison 
data for materials treated in vacuo and in gas 
atmospheres. Subjects discussed are ferrous 
base and high temperature alloys, refractory 
metals, nuclear power plant materials, elec- 
tronic materials, powder metallurgical prod- 
ucts, and vacuum coatings for metals. This 
section enables the reader to see what some of 
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the advantages of vacuum processing are and 
gives him an insight into its potentiatities. 

Three chapters are devoted to thermody- 
namics and kinetics. One of these is some- 
what elementary in covering very briefly such 
things as the Phase Rule, LeChatlier’s Princi- 
ple, the Gas Laws of Thermodynamics, etc. 
Also in another of these three chapters the 
figures are poorly reproduced. In general, 
however, this part of the text is useful. 

The chapter devoted to electron bombard- 
ment melting techniques is somewhat disap- 
pointing. It is written in such a general way 
that it is of very little practical value except 
to the person almost completely unfamiliar 
with the technique. A minor oversight is the 
lack of captions for some of the figures. 

Not all of the chapters contain references; 
however, the editor has compiled a rather com- 
plete and useful bibliography to supplement 
the text references. There are author and 
subject indexes. 

L. SMitH 
The Franklin Institute Laboratories 


ANALYSIS AND CONTROL OF NONLINEAR Sys- 
TEMS, by Y. H. Ku. 360 pages, diagrams, 
6X9 in. New York, The Ronald Press 
Co., 1958. Price, $10.00. 


This book is an outgrowth of a graduate 
course in nonlinear circuit analysis given at the 
Moore School of Electrical Engineering of the 
University of Pennsylvania by Professor Ku 

With the dawn of automation thundering 
all around us it becomes difficult to make a 
wise decision as to which books will be sig- 
nificant additions to one’s technical library 
and not just some more clutter. Professor 
Ku’s book is such a significant addition. 

Ever since the installation in the early fifties 
of Typhoon, the tremendous analog computer 
at the Naval Air Development Center, and 
smaller analog computers elsewhere, the phase- 
plane analysis has been tremendously useful 
in presenting analog or actual equipment 
results. 

In Professor Ku’s book a new approach has 
been developed based on the author’s method 
of space trajectory in m dimensions which ex- 
tends the phase-plane method to phase space. 
This is essentially a multidimensional plct al- 
lowing solution of higher order nonlinear dif- 
ferential equations than the second. 
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For the first time we have a book which 
provides a general and powerful method of 
studying the behavior of physical systems with 
many degrees of freedom. The method of 
simultaneous phase-plane equations demon- 
strates that a physical system with nonline- 
arities anywhere and with any degree of com- 
plication can be treated in a comprehensive 
manner. This method is a basic approach 
which extends nonlinear analysis to feedback 
systems with nonlinearities. Chapter 11 on 
analysis of nonlinear servomechanisms and 
Chapter 12 on nonlinear control of feedback 
systems are devoted to this new field of study. 

The preceding ten chapters carefully pre- 
pare the reader for the trenchant final chap- 
ters. Both analytical and topological methods 
of solving nonlinear differential equations are 
presented. Physical examples in both me- 
chanical and electrical systems which involve 
nonlinearities are worked out in detail. There 
is a particularly intriguing example in nuclear 
reactor kinetics at the end of Chapter 9 which 
is an analysis of high order nonlinear systems. 
Chapter 10 presents an analysis of multi-loop 
nonlinear systems and discusses the method of 
simultaneous phase-plane equations. This is 
the last stage of preparation for the final chap- 
ters. To aid the reader in acquiring technical 
dexterity ten well chosen examples are given 
at the end of each chapter for the reader to 
solve. 

After a leisurely and delightful journey 
through the ten chapters (delightful because 
of the beauty and simplicity of presentation), 
the reader is ready to partake of the piéce de 
résistance of the final chapters. In the chap- 
ter on analysis of nonlinear servomechanisms 
Professor Ku discusses responses in linear ser- 
vomechanisms, nonlinear elements for improv- 
ing servo performance, contactor servomecha- 
nisms, dual mode operation of nonlinear servo- 
mechanisms, multiple mode operation in a 
high order system, stability criterion for con- 
tactor servomechanisms, back lash in servo- 
mechanisms and nonlinear control systems for 
wide range input signals. 

In the final chapter he discusses the theory 
of nonlinear control of feedback systems, sys- 
tems with nonlinear control in the feedback 
branch, nonlinear feedback control systems 
with third and higher order transfer functions, 
systems with nonlinear control in the forward 
branch, nonlinear forward-control system with 
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integral feedback, nonlinear forward-control 
system with derivative feedback, feedback 
system with two nonlinear control functions, 
feedback system with two linear networks and 
two nonlinear control networks, nonlinear feed- 
back control to improve transient response of 
a system, nonlinear control in a torpedo sys- 
tem, stability study of higher order nonlinear 
control systems and applications of digital and 
analog computers in the design of nonlinear 
control systems. 

In addition to writing an invaluable text, 
Professor Ku has compiled a remarkably com- 
plete bibliography extending from 1860 to 
1958. This bibliography completes the book. 

Epwarp A. KNOBELAUCH 
Burroughs Corporation Research Center 


So_ip STATE Puysics, VOLUME 7, edited by 
F. Seitz and D. Turnbull. 525 pages, dia- 
grams, 6 X 9 in. New York, Academic 
Press Inc., 1958. Price, $12.00. 


The seventh volume of the SSP series main- 
tains the high standards of its predecessors, in 
the choice of topics and the competence of the 
treatments. We review successively the vari- 
ous subjects dealt with in this volume. 

1. Thermal Conductivity and Lattice Vi- 

brations Modes (pp. 1-99) P. G. Clemens 

After a brief introduction to this impor- 
tant subject, the author devotes about a third 
of his exposé to an excellent presentation of 
the necessary background, namely the theory 
of quantized lattice waves, normal modes, the 
various types of interactions and to the ap- 
proximate solution of Boltzmann's transport 
equation. He then treats successively the 
thermal conductivity of (a) Non-Metals (at 
high and low temperatures, imperfections, size 
effect etc.) and of (b) Metals, Alloys and Semi- 
conductors (lattice and electronic components 
and their separation etc. . . .). This presen- 
tation exhibits the same qualities of clarity 
and order, shown by the author in his earlier 
monograph on a similar subject (on Low Tem- 
perature Thermal Conductivity, Handbuch 
der Physik, Vol. 14, Berlin 1956). 

2. Electron Energy Bands in Solids (pp. 

100-209) J. Callaway 

This is a central and difficult problem of 
Solid State Physics, and it will remain so for 
many years tocome. J. Reitz in SSP Vol. 1 
and T. O. Woodruff in Vol. 4 have already 
dealt with the fundamentals of the subjects. 


: 

2 

ES 

: 

3 

the 

: 

i 

ar 


256 


J. Callaway after a fifteen page critical sum- 
mary of the general theory, discusses succes- 
sively the alkali metals, the metals of groups 
II and III, the group 1V elements (and related 
semiconductors), the transition metals and 
the 4 other substances. This monograph is 
remarkable by the wealth and scope of the 
information provided, and is a must for any- 
one interested in the energy band theory of 
solids. 

3. The Elastic Constants of Crystals (pp. 
214-352) H. B. Huntington 

The classical theory of elastic constants 

is introduced along with the atomistic theories. 
An exposé of the various experimental methods 
is followed by an excellent and most useful set 
of tables of measured elastic constants, of vari- 
ous metal groups, rare gas solids, ferromagnetic 
materials etc. The last part is a lucid and 
thorough discussion of the dependence of 
elastic constants on Temperature, Pressure, 
Composition Phase Changes and Relaxation 
Phenomena. 

4. Wave Packets and Transport of Elec- 
trons in Metals (pp. 353-378) H. W. 
Lewis 

After a brief restatement of the current 

theory of electronic transport, the author ana- 
lyzes some of its shortcomings and their con- 
sequences. Emphasis is put on the unadia- 
batic nature of the scattering processes in 
metals. 

5. Study of Surfaces by Using New Tools 
(pp. 379-425) J. A. Becker 

With the growing importance of the 

physics of surface phenomena, Becker's pres- 
entation of three new and powerful experi- 
mental tools of surface investigation will be 
welcomed by many physicists. The three 
tools are the ion gauge, the FEM (field emis- 
sion microscope) and the FIM (field ion micro- 
scope). The ion gauge measures the rates of 
adsorption and desorption, the FEM yields 
the dependence of the work function on the 
crystallographic planes and the FIM allows 
the location of individual atoms on metal 
surfaces. 

6. The Structures of Crystals (pp. 426-503) 
A. F. Wells 

This is an illuminating and up-to-date 

survey. The author is to be commended for 
the clarity of his exposé and especially for the 
excellent diagrams and tables. It is certain 
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that this article will be frequently and profita- 
bly consulted by many physicists and chemists. 
HENRI AMAR 

The Franklin Institute Laboratories 


ROCKET PROPELLANTS, by Francis A. Warren. 
218 pages, illustrations, 6 X 9 in. New 
York, Reinhold Publishing Corp., 1958. 
Price, $6.50. 


With the coming of age of astronautics many 
books have appeared on satellites and space 
flight. In most of these some bits of informa- 
tion have been given concerning propulsion 
systems which permit the payloads to go into 
orbit or beyond the Earth’s gravitational field. 
In these books we may find reterences to liquid 
and solid propellants, to rocket motors, to 
pumps, to igniters but always the treatment is 
meager almost to the point of frustration. 

Now appears Rocket Propellants written by 
Francis A. Warren who, in working for the 
Southwest Research Institute, coordinates re- 
search on propellants and their combustion. 
From an excellent background in propellants 
he has put together a fascinating story on why 
rockets work. The author claims this is writ- 
ten‘. . . to provide technical men with basic 
information on the materials being used to 
propel the rockets and missiles of today, ....” 
The author has capably fulfilled his commit- 
ment not only to the technical personnel but 
also to those non-technical people who find 
propellants absorbing and would like to know 
about them. 

In the literature we find many references to 
the various types of propellants, the grain 
shapes typical of the state of the art today, the 
ingredients of both the solid and liquid propel- 
lants, the processing the “‘pros’’ use in produc- 
ing these fuels and the performance character- 
istics of propellants. Here for the first time 
we find this material collated in a book. 

When we watch a big rocket being launched 
and see the plume of fire roaring down on the 
deflector little does the layman realize the tre- 
mendous volume of research which has been 
successfully carried through to make possible 
the precise burning of the propellants. Chap- 
ter V, The Burning of Propellants, is a classic 
treatment of what happens when the rocket is 
ignited. Temperature, burning rate, velocity 
of exhaust are some of the topics discussed in 
this chapter. The author in comparatively 
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few pages reveals much of the burning process. of safety in manufacturing and handling pro- 
Ignition and igniters are given a chapter. _ pellants and finally there is a projection into 
This is one rarely mentioned in conventional — the future. 
books yet this reviewer finds it one of the most The author has put together a readable and 
fascinating. Extreme ingenuity has been ex- useful book which will be much in demand in 
ercised by designers to insure the ignition of _ the engineer's library. This book was long 
the rocket. This chapter shows what must overdue and it may well set the pattern for 
be done to realize the burning state. others to follow. 
The author concludes with chapters on solid 
and liquid propellants, the important element 


I. M. Levitt 
The Franklin Institute, Planetarium 
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Price, $2.55 (paper). 

E.uiptic INTEGRALS, by Harris Hancock. Unabridged and unaltered republication of the 
first edition. 104 pages, diagrams, 54 X 84 in. New York, Dover Publications, Inc., 
1959. Price, $1.25 (paper). 
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THE FOUNDATIONS OF EUCLIDEAN GEOMETRY, by Henry George Forder. Unabridged and 
unaltered republication of the first edition. 349 pages, diagrams, 54 K 84in. New York, 
Dover Publications, Inc., 1959. Price, $2.00 (paper). 


THE THEORY OF THE POTENTIAL, by William Duncan MacMillan. Unabridged and unaltered 
republication of the first edition. 469 pages, diagrams, 5} X 84 in. New York, Dover 
Publications, Inc., 1959. Price, $2.25 (paper). 


AN INTRODUCTION TO THE GEOMETRY OF N DiMENSIONS, by D. M. Y. Sommerville. Un- 
abridged and unaltered republication of the first edition. 196 pages, diagrams, 54 X 84 in. 
New York, Dover Publications, Inc., 1959. Price, $1.50 (paper). 


THE GENETICAL THEORY OF NATURAL SELECTION, by Ronald A. Fisher. Second revised 
edition of the work first published in 1929. 291 pages, diagrams, 5 X 8in. New York, 
Dover Publications, Inc., 1959. Price, $1.85 (paper). 

Evasticity, PLasticiry & STRUCTURE OF MATTER, by R. Houwink. Unabridged republica- 
tion of the second edition. 368 pages, diagrams, 5$ X 8in. New York, Dover Publica- 
tions, Inc., 1959. Price, $2.45 (paper). 

THE THEORY OF THE PROPERTIES OF METALS AND ALLOoys, by N. F. Mott and H. Jones. 
Reprinted from corrected sheets of the first edition. 326 pages, diagrams, 54 X 84 in. 
New York, Dover Publications, Inc., 1959. Price, $1.85 (paper). 

STUDIES ON THE STRUCTURE AND DEVELOPMENT OF VERTEBRATES, by Edwin S. Goodrich, 


Two-volume unabridged and unaltered republication of the one-volume first edition. 837 
pages, diagrams, 5} K 8 in. New York, Dover Publications, Inc., 1959. Price for each 


volume, $2.50 (paper). 


GUIDE TO THE LITERATURE OF MATHEMATICS AND Puysics INCLUDING RELATED WoRKS ON 
ENGINEERING SCIENCE, by Nathan Grier Parke III. Revised and expanded version of 
the first edition. 436 pages, 5} X 8 in. New York, Dover Publications, Inc., 1958. 


Price, $2.49, (paper). 
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New Powder Pattern Technique. — 
A new method for delineating the do- 
main structure at the surfaces of ferro- 
electric crystals was described recently 
in a paper presented to the American 
Physical Society by G. L. Pearson and 
W. L. Feldmann of Bell Telephone 
Laboratories. It has been used to de- 
lineate domain structures in great de- 
tail on a wide variety of crystals in 
which domains have never been ob- 
served before. 

According to its discoverers, the 
method should be applicable to the 
study of piezoelectric and pyroelectric 
charged surfaces as well as ferroelec- 
tric materials. 

The new technique uses colloidal sus- 
pensions of electrostatically charged 
powders in an insulating organic liquid. 
When a few drops of this suspension 
are applied to the face of the crystal, 
the charged powder is immediately at- 
tracted to the ferroelectric domains 
carrying an opposite charge, and covers 
their entire area. 

The most effective materials for this 
work are commercial spray-grade sul- 
fur and red lead oxide, each suspended 
in hexane. The sulfur deposits on 
negatively charged domains, while the 
lead oxide deposits on the positively 
charged domains. Each of the sus- 
pensions is applied separately ; the sec- 
ond is not applied until the hexane in 
the first has evaporated. If the two 
suspensions mix on the crystal sur- 
face, each may lose itscharge. In this 
event, no pattern is formed. The 
powders are fixed in place indefinitely 
by their electrostatic charge after the 
hexane evaporates. 

The yellow sulfur and red lead oxide 
provide brightly colored delineation of 
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the positive and negative domains, and 
the pattern shows great detail. 

A dispersion of a cross linked poly- 
mer derived from polystyrene can also 
be used as the negatively charged col- 
loid, in place of the lead oxide. This 
can be dyed any desired color with an 
oil soluble dye. 

The basis for the success of the new 
technique lies in the fact that although 
the colloid as a whole is electrically 
neutral, individual particles acquire 
a diffuse, double-layer charge when 
brought in contact with the liquid. 
Under the influence of the ‘“built-in’’ 
electric field, the colloidal particles are 
attracted either to the positive or nega- 
tive domains depending on the orien- 
tation of their dipolelayers. For these 
experiments, an insulating liquid of low 
viscosity and low dielectric constant 
such as hexane is desirable, so that 
the charged particles are free to move 
toward the ferroelectric domains under 
maximum electrostatic attraction. 

Powder patterns of an entirely differ- 
ent nature result if a colloid of finely di- 
vided barium titanate in hexane is 
used instead of sulfur and lead oxide. 
The BaTiO; particles carry no net elec- 
tric charge, but they do have a very 
high dielectric constant. Thus, when 
they are placed on the surface of a fer- 
roelectric crystal containing positive 
and negative domains, the particles 
polarize in the external electrostatic 
fields and are deposited at domain 
boundaries where large field gradients 
exist. The resulting powder patterns 
outline the domain boundaries well, 
but do not differentiate between posi- 
tively and negatively charged domains. 

This new technique has provided the 
first information available on the do- 
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main structure of a number of crystal- 
line ferroelectric materials, including 
triglycine sulfate and guanidinium gal- 
lium selenate hexahydrate. Previous 
optical determinations which have 
been made on Rochelle salt, and do- 
main etching procedures on barium ti- 
tanate have been confirmed by the 
new method also. 

Thus, a simple, easily employed 
technique joins the group of tools 
available to the physical chemist and 
physicist in their continuing search for 
knowledge of the ultimate structure 
of matter. 


The Proceedings of the Royal Insti- 
tution of Great Britain.—The Proceed- 
ings of the Royal Institution are now 
being published three times a year in- 
stead of annually. The Proceedings 
contain accounts of the well-known 
Friday Evening Discourses and one ad- 
vantage of the new procedure will be 
that the Discourses will appear in print 
soon after their delivery. 

The Evening Discourses at the 
Royal Institution were started by 
Michael Faraday in 1826 and for over 
130 years they have maintained the 
highest standards of popular exposi- 
tion. The Discourses are famous, not 
only because they are given by leading 
men of science, but also because it is 
traditional to illustrate the subject 
by experiments and demonstrations 
whenever possible. 

The Proceedings therefore contain 
authoritative accounts of research and 
learning, written for the non-specialist 
and covering many branches of science 
and other subjects. In these days of 
specialization, such accounts enable 
the reader to get a broad survey of the 
most important advances in a wide 


field. 


Moisture Balance with Built-in Au- 
totransformer.——-A new and improved 
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moisture balance for fast and highly 
accurate determinations of moisture 
content in a wide variety of materials 
has been designed by Central Scientific 
Company, Chicago. Principal feature 
of the new model is its built-in auto- 
transformer which regulates voltage 
automatically and provides more con- 
venient temperature control than pre- 
vious moisture balances in which the 
autotransformer was a separate unit. 

Among the materials for which the 
balance may be used with a high de- 
gree of accuracy and reliability are: 
nonfat milk powder, dry milk, soy- 
bean products, pharmaceuticals, chem- 
icals, corn products, leather, fiber, la- 
tex, paper pulp, soap and detergents, 
food products, ceramic powder, clays, 
resins, paint pigments, soil, raw meat, 
fertilizer, dextran, yeast, plastics, lac- 
tose, plant seeds, and cotton. Prod- 
ucts that do not change their chemical 
structure or decompose while losing 
water are most suitable for testing. 

Simple to operate by even an inex- 
perienced laboratory technician, the 
moisture balance is capable of reducing 
from one-half hour, or more, to only a 
few minutes the time required for ac- 
curately testing the moisture content 
of a material. 

According to Central Scientific, the 
balance performs both drying and 
weighing operations simultaneously, 
thus eliminating time-consuming weigh- 
ing of samples on an analytical bal- 
ance, with subsequent drying in a 
vacuum oven and then reweighing. 

To operate the Cenco moisture bal- 
ance, a 5-gram sample (exact weight 
need not be known) is placed in the 
balance pan and the indicator crank is 
set for 100 per cent moisture. The 
hood is then lowered, the lamp switch 
turned on, and the infrared heat quickly 
penetrates deeply into the sample. A 
few minutes later, with moisture driven 
off, the percent of moisture lost is read 
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directly on an easy-to-read calibrated 
scale. 

The improved unit also comes 
equipped with disposable weighing 
dishes, made deeper than former mod- 
els. These provide a better flexibility 
of loads and enable the operator to use 
larger liquid samples, or bulkier dry 
samples. The dishes are ridged to 
keep liquid samples uniformly distrib- 
uted in the weighing pan and elimi- 
nate spillage. 

The built-in magnetic damper re- 
duces pan vibrations and speeds weigh- 
ings. An extra torsion wire for 25 
gram samples is supplied for handling 
heavier loads. 

The new Cenco moisture balance is 
housed in a completely redesigned and 
compact case of solid cast aluminum 
with Cenco blue wrinkle enamel. It 
can be used wherever a 115-volt, 60 
cycle outlet is available. 


In-Circuit Transistor Tester.—The 
first transistor tester capable of check- 
ing the performance of transistors 


while they are connected within their 
circuits, has been developed by Philco 
Corporation, and a commercial model 
will be available in the near future. 
In-circuit transistor testing is ex- 
tremely important since the method 
saves trouble shooting and servicing 
time without turning on power in the 
equipment. 

Other commercially available tran- 
sistor testers cannot perform in-circuit 
transistor measurements since they are 
unable to differentiate between normal 
input and output signals generated 
by the transistor and spurious signals 
appearing at either the input or out- 
put terminals. These spurious signals 
arise from the sneak paths propagated 
by the circuitry surrounding the tran- 
sistor under test. 

Low impedance methods are em- 
ployed in both the input and output 
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circuits to nullify the effects of the ex- 
ternal circuitry. 


Cathode Follower Probe Improves 
Accelerometer Performance.—A new 
type of Cathode Follower Probe, re- 
quiring no separate power, has been 
developed by Columbia Research Lab- 
oratories. Composed entirely of pas- 
sive circuitry, this new device provides 
an improved solution to the problem 
of using long cables up to 100 feet 
or more without loss of accelerometer 
sensitivity. 

Basically, the device functions as 
the first stage of the Cathode Follower 
and because of its small size and light 
weight (? X 4-in. diameter and 3 oz.) 
the unit can be mounted close to the 
accelerometer. This has the effect of 
bringing the Cathode Follower out to 
the transducer. In low g-level studies 
this affords the possibility of using 
very short low-noise cable going into 
the probe, thereby providing a means 
of increasing sensitivity as much as 
four times with some accelerometers. 

The Probe connects directly to the 
input jack of the Cathode Follower 
and requires no circuit adjustment to 
install. It is detachable, extremely 
small and rugged, free from micro- 
phonics and capable of withstanding 
temperatures up to 600° F. 


Acoustoelectric Effect in Germa- 
nium.—aA direct current can be in- 
duced in a semiconductor under cer- 
tain conditions by passing acoustic 
waves through it. This effect, which 
has been studied extensively by Gab- 
riel Weinreich of Bell Telephone Labo- 
ratories, is known as the “‘acoustoelec- 
tric effect.” 

An acoustic wave traveling through 
a semiconductor tends to produce a 
net force acting on the charge carriers 
in the direction of propagation of the 
wave. To offset this force, an electric 
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field is induced in the semiconductor 
material, and this field can be meas- 
ured by attaching contacts along the 
length of the sample. 

Since this effect depends on some 
bunching of the carriers, it is normally 
reduced greatly by their mutual elec- 
trostatic repulsion. There are some 
situations, however, in which more than 
one type of carrier is present, and it is 
possible to bunch each type separately 
without developing a net space charge. 
Such is the case with n-type germa- 
nium, in which the electrons belong 
to a band which has a number of 
distinct energy minima or “‘valleys.”’ 
The electrons in the vicinity of these 
various valleys react differently to the 
acoustic wave, giving rise to the acous- 
toelectric effect. The magnitude of 
this effect gives a measure of the inter- 
valley scattering time, that is, the 
mean transition time between the sepa- 
rate groups of states in the conduction 
bands which lie near the minima of 
energy (valleys). 

Experimental verification of the 
acoustoelectric effect has been ob- 
tained by passing acoustic waves at 
frequencies of 20 mc. and 60 me. 
through a bar of n-type germanium, 
and measuring the voltage induced 
along the bar. To minimize stray 
effects, the acoustic waves are modu- 
lated with an audio frequency, thus 
producing an audio frequency induced 
voltage rather than a d-c. voltage. 

A study of intervalley scattering 
time has led to a better understanding 
of semiconductor materials. It has 
been found that this scattering is due in 
part to phonon scattering and in part 
to impurity scattering. These two 
effects may be separated, as the pho- 
non effect becomes greater at higher 
temperatures (above 100° K.) while 
the impurity effect is predominant at 
temperatures below about 60° K. Be- 
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tween 60° K. and 100° K., the two 
effects are mixed. 


Ferroelectric Polarization Reversal 
in Barium Titanate.—Recent experi- 
ments on single crystal barium titanate 
have demonstrated for the first time 
that complete ferroelectric polariza- 
tion reversal can be accomplished 
under an applied electric field, with a 
single expanding domain which grows 
through extensive sidewise motion of 
the 180° domain walls. These experi- 
ments were described in a paper by 
Robert C. Miller of Bell Telephone 
Laboratories, which was presented to 
the Conference on the Physics of Di- 
electrics at Moscow. 

Very little work has been reported 
previously on sidewise motion of 180° 
domain-walls in ferroelectrics. These 
180° domain-walls are the boundaries 
between two adjacent areas or do- 
mains of opposite polarization. Side- 
wise motion involves the movement of 
boundaries in a direction perpendicu- 
lar to the dipole moment of the do- 
main. In Dr. Miller’s work, a single 
crystal of BaTiO; is mounted with 
liquid electrodes applied to either side, 
and is subjected to an external applied 
field. The electrodes used were silver 
or platinum rods 1.0 mm. in diameter, 
with a small hole drilled in one end. 
An aqueous lithium chloride solution, 
saturated at room temperature, was 
placed in the hole as the electrolyte. 

Partially switched crystals were 
studied by removing them from the 
apparatus and etching them fora short 
time in weak hydrofluoric acid. The 
positive dipole ends etch preferen- 
tially ; examination with a microscope 
after etching allows one to observe the 
domains visually. 

With the liquid electroded crystals, 
the number and size of the reversed 
domains in the partially switched crys- 
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tals was of great interest. In numer- 
ous cases, a single reversed domain 
comparable in size to the total elec- 
trode area was observed. Many of 
the domains appeared to nucleate at 
obvious crystalline imperfections. 

In a further refinement of the tech- 
nique, Miller took advantage of this 
nucleation phenomonon by “‘manufac- 
turing’ a nucleation site. This site 
was a sandblasted dimple about one- 
third the crystal thickness deep, and 
with a diameter about one-fourth that 
of the electrode. 

When a dimpled sample is placed in 
the liquid electrode holder, nucleation 
of reversed domains is most probable 
in the dimpled region because it is sub- 
jected to a high field compared to that 
of the surrounding area. Once a nu- 
cleation has occurred in the dimpled 
region, the domain grows out from the 
dimple through sidewise wall motion. 
The single domain is roughly square, 
as has been shown by successive etch- 
ing of the partially switched crystal. 
In this way, complete reversal of polar- 
ization has been accomplished for the 
first time with a single expanding do- 
main, starting from a predetermined 
point on the crystal. 

There is some evidence to indicate 
that impurities have an effect on the 
velocity of the domain wall motion, 
but the evidence is meager. The ex- 
periments performed here indicate that 
the velocity is given by the equation: 


for over three decades of velocity. 
The velocity for the field E= & is v,, 
and is of the order of 10 cm/sec. The 
activation field varies from 2000 to 
4000 v/cm. 

It appears that under these condi- 
tions, this polarization reversal is anal- 
ogous to magnetization reversal in 
some ferromagnetic materials. 
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Television System for High Altitude 
Solar Observation.—<A special televi- 
sion system that will permit astrono- 
mers on the ground to aim and focus a 
telescope suspended from a_ balloon 
fifteen miles above the Earth is being 
developed by scientists of the Radio 
Corporation of America. The tele- 
vision equipment, comprising a spe- 
cially designed slow-scanning airborne 
camera and transmitter as well as a 
ground monitoring system, is to be 
used in high-altitude solar observation 
by a Princeton University team under 
Dr. Martin Schwarzschild of the De- 
partment of Astronomy. 

In the experiments, scheduled for 
next summer from a base in Min- 
nesota, a 12-in. reflecting telescope, 
known as the ‘‘Stratoscope,”’ will be 
carried to altitudes of 80,000 ft. for 
observation of the sun’s surface. At 
these altitudes, observations may be 
made with little or no interference from 
the Earth’s atomosphere, whose den- 
sity and turbulence limit clear viewing 
from ground observatories. 

The special RCA television camera 
will be mounted on the “Stratoscope”’ 
and linked optically with a motion 
picture camera that obtains a perma- 
nent photographic record of the im- 
ages appearing in the 12-in. reflector. 
Through the airborne transmitter, also 
mounted on the telescope structure, 
the television camera will transmit to 
the ground a second-by-second view of 
what the motion picture camera sees. 
Watching the results on television 
monitor screens, the astronomers at 
the ground base will control the aim 
and focus of the telescope by means of 
a remote control system that also is 
being developed for the project. 

The ‘‘Stratoscope” program, initi- 
ated and directed by Dr. Schwarzs- 
child, produced a number of clear and 
detailed views of the solar surface dur- 
ing two unmanned flights in 1957. At 
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that time, however, the system did not 
include television and remote control 
equipment to keep the telescope in 
proper focus and aimed directly at the 
sun at all times. As a result, only a 
portion of the photographic record 
made by the motion picture camera 
provided the desired information. Ac- 
cording to Dr. Schwarzschild, the de- 
velopment of the new television and 
control systems should result in full 
control of the “Stratoscope’’ through- 
out its flight, producing a complete 
photographic record of useful informa- 
tion for the astronomers. 

The special RCA television system, 
being developed under the supervision 
of Dr. V. K. Zworykin, Honorary Vice 
President of RCA, is designed to pick 
up and transmit to the ground one 
picture per second up to a range of 
150 miles. One picture per second will 
provide all of the information needed 
by the astronomers for control pur- 
poses, and will at the same time permit 
the use of low power equipment trans- 
mitting over a far narrower bandwidth 
than that required for conventional 
television broadcasting. 

The special television camera, which 
will be mounted beside the motion pic- 
ture camera on the telescope structure, 
employs a 1-in. RCA Vidicon pick- 
up tube and transistorized circuits. 
Measuring only 8 X 6 X 5 in., it will 
be surrounded by a pressurized case for 
high altitude operation. Thecompact 
frequency-modulation (FM) transmit- 
ter, together with its battery power 
supply, also will be contained in a 
pressurized case. <A third element in 
the airborne system isa special antenna 
arrangement contained in a shock-re- 
sistant foam package suspended below 
the telescope. The antennas willsend 
television and telemetering signals to 
the ground station and will receive 
from the ground the control signals for 
the focusing and aiming devices. 
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The television camera will be linked 
optically to the motion picture camera 
by a mirror system. In operation, a 
shutter arrangement will cause the 
telescope image to be reflected briefly 
once per second to the television cam- 
era, which will transmit to the ground 
during the remainder of the second a 
picture identical to that recorded at 
the same instant by the film camera. 

The ground equipment will include 
a television receiver and three moni- 
tors, one of which will be used for 
photographing the transmitted images 
for later checking against the perma- 
nent film record. The entire televi- 
sion system will operate at a frequency 
of about 225 megacycles, in the ultra- 
high-frequency (UHF) band, on a 
bandwidth of 200 kilocycles. 


New Design for School Labora- 
tories.—Fifteen seconds after the flip 
of aswitch, a physics classroom labora- 
tory was converted to a fully-equipped 
biology laboratory in a recent demon- 
stration at the Chicago Undergraduate 
Division of the University of Illinois. 
Physics Professors Ora L. Railsback 
and Harold M. Skadeland presented 
the revolutionary design in school lab- 
oratories in the form of a full-scale, 
working model at Navy Pier. The 
key to the quick change is a built-in, 
automatically rotated turntable at 
each student station which delivers 
one set of lab equipment “‘on stage” 
while a previous set disappears into 
the wall. All units can be rotated at 
once or each one singly from a centrally 
operated electronic panel. 

The new design thus allows a single 
laboratory to serve a variety of study 
courses——physics, biology, geology and 
others. This is a welcome economy 
for schools and colleges facing inflated 
building costs, Railsback and Skade- 
land believe. 

Its most obvious advantage is that 
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while one section is in use in the labo- 
ratory, the other sections can be loaded 
with equipment from a_behind-the- 
scenes storage room. ‘This contrasts 
with conventional school laboratories 
that often are taken out of use for a 
full day each week to be re-equipped. 

The professors presented a small 
scale model for the first time last June 
to the building committee of the Amer- 
ican Physics Teachers’ Association. 
Committee members asked tough ques- 
tions about the new design, but were 
satisfied with the answers. The po- 
tentials for saving time, space and 
money were obvious. 

Professors Railsback and Skadeland 
then went ahead with the full-size, work- 
ing model. Built by Hamilton Manu- 
facturing Company of Two Rivers, 
Wisconsin, the new lab, described as 
“still experimental’ by its designers, 
features two distinct types of student 
stations and turntables, each with its 
particular set of advantages. 

One grouping consists of a single, 
long laboratory bench with the series 
of turntables built into the wall above 
it. This permits a massive loading 
and unloading of units behind the wall. 

In the other grouping the units are 
set up individually at right angles to 
the wall. The turntables are loaded 
from within the classroom itself, not 
from a hidden storage area. The time 
factor in loading was not seen as criti- 
cal in the use of these units where the 
labs may be vacant for longer periods. 
They are particularly adaptable to re- 
search work of a continuing nature. 

The relative merits of both types 
of units are yet to be determined 
through further study, according to 
Dr. Railsback. 

With either type unit, equipment 
used in experiments of a continuing 
nature does not have to be dismantled 
at the close of each laboratory period 
and reassembled at the start of the 
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next. Instead, the equipment simply 
disappears into the turntable between 
sessions. 

It is conceivable that three continu- 
ing experiments could be carried on at 
each student station by three different 
student classes. 

The new design also remedies the 
cluttered, disorganized appearance, usu- 
ally found in advanced and gradu- 
ate class labs. And, of course, the 
speed with which the lab can be con- 
verted provides students and teachers 
alike with more total time for ex- 
perimentation. 

Each station is equipped with water 
faucet and drain and electrical outlets 
with variable voltage outputs. The 
elaborate electrical circuits feed into a 
large, futuristic ‘“‘switchboard.”’ 

“Rising student enrollments are go- 
ing to tax the ingenuity of every scho- 
lastic institution to make the maxi- 
mum use of available space,’ Dr. 
Railsback declared. ‘‘We think our 
new laboratory design is a positive 
step in that direction.” 


New TV Lens.—A new television 
lens—technically known as a Studio 
Varotal Zoom Lens—has been devel- 
oped in Britain as the result of close 
collaboration between the BBC and 
Taylor, Taylor and Hobson (a division 
of Rank Precision Industries Ltd.), 
pioneers in this field. 

It enables cameramen to make grad- 
ual changes between panoramic views 
and magnified close-ups without mov- 
ing the camera and without losing any 
quality of definition. 

Designed specifically to cater to the 
peculiarities of the television camera, 
and said to be a great technical ad- 
vance, the lens was first used to tele- 
vise Queen Elizabeth’s Christmas Day 
address. 

The full value of zoom lenses can 
only be obtained if they provide as 
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good quality definition as normal lenses 
which they replace. 

This requirement presents severe 
optical and mechanical difficulties. 
Taylor Hobson overcame these diffi- 
culties for outside television broad- 
casting some years ago. They pro- 
duced a zoom lens—the Outside Broad- 
cast Varotal model—that in the past 
six years has become internationally 
accepted. 

When similar results are required in 
the restricted space of an indoor tele- 
vision studio, the problems increase 
greatly. An Outside Varotal is not 
generally suitable for these conditions. 
The new Studio Varotal was designed 
to meet the worldwide demand for a 
zoom lens for studio use. 

This is suitable for use with both 
image orthicon and vidicon camera 
tubes, interchangeable rear units being 
provided to suit the two formats. The 
focal length range of the Studio Varo- 
tal is 2} in. to 8 in. when fitted on 
image orthicon cameras and 2.25 cms. 
to 8 cms. on vidicon cameras. 

The lens will focus down to 5 ft., 
which is unusually close for a zoom 
lens. These features have been pro- 
vided while maintaining adequate rel- 
tive aperture and high standard of 
definition. 

The optical design of the Studio 
Varotal has been specially developed 
to yield the type of performance de- 
manded by television transimssion 
channels. The correction of all lens 
aberrations matches the spectral sensi- 
tivity of the television camera tube, 
and aberrational compensation has 
been provided for errors introduced by 
the face of the tube. 

From the mechanical point of view, 
the lens has been designed to fit exist- 
ing camera turrets with a minimum 
of adaptation. The three controls— 
focus, zoom and iris—take the form of 
gear rings mounted towards the rear 
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of the lens, facilitating the fitting of 
either manual or servo drives. 

Camera controls for the lens used 
during the Queen’s broadcast were 
specially developed by the Engineer- 
ing Division of the BBC. The lens is 
operated by a capstan on the right- 
hand side of the camera. This rotates 
for focusing and slides backwards or 
forwards for zooming, enabling the 
cameraman to perform the whole oper- 
ation with one hand. 


Oscar Fisher Company Unveil New 
Automatic X-Ray Processor.—-The 
Oscar Fisher Company, of Newburgh, 
N. Y., have developed a new auto- 
matic processor, the Fisher FaX-Ray 
Processor. 

The new unit, entirely automatic in 
operation, develops, fixes, washes and 
dries X-ray film of all sizes, up to and 
including 14 in. by 17 in. One of the 
unique features of this new apparatus 
is that the film transport is accom- 
plished without making use of cum- 
bersome, old-fashioned hangers. 

Present methods of processing X-ray 
film require a minimum of one hour. 
The Fisher FaX-Ray Processor will 
make available to doctors completely 
dry X-rays in a matter of only 7 min- 
utes. In critical surgical operations 
and emergency treatments in the na- 
tion’s hospitals, the increased rate of 
processing, due to the use of a FaX- 
Ray Processor, will undoubtedly help 
the pathologist to save lives and re- 
lieve pain and suffering. 

‘The FaX-Ray Processor also reduces 
the dosage required for X-ray diagno- 
sis, allowing doctors ample time to ex- 
amine films at leisure rather than sub- 
jecting patients to continuous X-ray 
or fluorography. 

The FaX-Ray Processor is fabri- 
cated entirely of type 316L stainless 
steel. This is the only one of 42 differ- 
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ent types that satisfactorily resists the 
corrosive attacks of photographic 
chemicals. All seams Heliarc 
welded and surfaces passivated. 

X-ray films are developed in the 
FaX-Ray Processor by the Fisher 
method of agitation-immersion. 
Washing, too, isspeeded up by aspray- 
thru-air system, after which X-ray 
films are dried by turbulent, filtered 
air and delivered automatically in only 
7 minutes. 

The Fisher FaX-Ray Processor can 
be operated by anybody and does not 
require photographic or technical train- 
ing of any kind. Because the unit is 


entirely automatic, the operator after 
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ten minutes of instruction, has only to 
check the proper temperatures of the 
solutions before setting the FaX-Ray 
Processor into operation. The 12-film 
automatic feed magazine is loaded in 
a darkroom with exposed film, from 
where it is transported automatically 
through a pass-thru chute to the FaX- 
Ray Processor in an adjacent room. 
Placing the FaX-Ray unit outside the 
darkroom permits the operator to 
work at peak efficiency and accuracy 
under daylight conditions. 

Approximate dimensions of the 
Fisher FaX-Ray Processor are : length, 
6 ft.; height, 7ft.; width, 4 ft.; and 
weight about 500 Ib. (dry). 
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THE FRANKLIN INSTITUTE LABORATORIES 
FOR RESEARCH AND DEVELOPMENT 


In the fields of engineering and the physi- 
cal sciences, a competent and versatile 
staff of several hundred scientists and en- 
gineers, working with modern equipment 
in a creative climate, is trained to bring a 
fresh scientific approach to the solution of 


industrial problems. 


THE FRANKLIN INSTITUTE OF THE STATE OF PENNSYLVANIA 
BENJAMIN FRANKLIN PARKWAY AT 20TH ST. PHILADELPHIA 3, PA. 
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